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1.0 INTRODUCTION 


Noise generated bv fan and compressor stages is a complex spectrum of 
tone and broadband sources produced by several source mechanisms. One major 
source of fan noise is the interaction of the nonuniform, unsteady rotor exit 
flow field with the downstream stator vanes. Rotor exit flow nonuniformities 
and unsteadiness produce fluctuations in flow velocity and flow angle into the 
stator, thereby generating unsteady fluctuating lift forces on the stator 
vanes. These rotor exit flow nonuniformities and unsteadinesses are pro uceo 
by the viscous wakes shed from the rotor blade trailing edges and by the sec- 
ondary flow vortices formed In the rotor blade passages. 

The combined effects of blade wake and hub and tip vortex formations 
result in a rather complex rotor exit velocity field, seen by the downstream 
stator vanes as a time-varying, periodic gust excitation with a complex wave- 
form. The resulting vane row unsteady lift exhibits a spectrum containing 
tones at integral multiples of rotor blade passing frequency (BPF) , as well 
as some broadband response, thought to be associated with the random unsteadi- 
ness and turbulence produced in the blade wakes and end wall (hub and tip) 
regions of the rotor exit flow. The resultant noise spectrum is directly 
related to the vane unsteady lift spectrum through the coupling of the 
unsteady lift energy to acoustic modes in the fan duct. 

It is thus clear that the vane row unsteady lift response and the result- 
ing noise spectrum is directly linked to the character (spatial and time vari- 
ations) of the combined wake/vortex flow field of the rotor. In order to be 
able to adequately predict this rotor/stator interaction and noise spectrum, 
a realistic model of a rotor exit wake/vortex flow field is required. 

The major objective of the program summarized in this report was to estab- 
lish and verify a rotor wake and vortex model for specific application to fan 
and compressor rotor-stator interaction noise generation. The model could 
also be used to evaluate stator vane forced vibration excitations. A primary 
requirement for this model was that it be, in its final computer code form, 
compact and computationally rapid so that it could easily be adapted as a sub 
program or subroutine to a larger fan noise prediction computer program. 

The approach selected for development of the rotor wake and vortex model 
was to analyze existing rotor exit flow wake and vortex data and evolve an 
empirical or semiempirical correlation. This correlation was tc be guided 
by parallel analytical studies aimed at identifying the appropriate similarity 
parameters in order to cast the correlation into a "universal" or general form. 
The specific approach taken was to develop the empirical correlations of rotor 
wake and vortex properties (such as wake centerline velocity defect, wake semi- 
width, etc. ) in a rotor-fixed (rotating) reference frame, then mathematically 
transform these results into a stationary, stator-fixed reference frame. 


r » 


. f j Th f p T es<int: Program comprises two major tasks. The first task consists 
f developing a mathematical model and the associated computer program to pre- 
dict viscous wakes and secondary flow vortex velocity fields downstream of ? a 
rotating o’ade row. This rash includes a 1 iterator*' search and” collation 

Saaiam a frf ?" A ? Xi “ ing rotor for use in the model development. 

Semiempirical models of the rotor wake mean flow and turbulence 1-vensitv char- 
acteristics have been developed using the data compiled from the literature 
survey, and a parallel analytical model of the viscous wake has been developed 

model^ of S rotor t exi^° 8re % S ° f ^ senilen,pi rical model - A simplified analvtical 
a?so been? secondary vortex formation and downstream development has 

aiso been incorporated into the rotor wake and vortex model. The resulting 

linked to thl° r r h° the staCor vane upwash velocity spectrum, and is 

linked to the steady aerodynamic characteristics of the fan rotor. 

In TalkV^h™ li** COnSls f s of evaluating the prediction procedure developed 
n Task I, through parametric studies to assess the behavior of the model when 

r ^ are r rled 3nd thr ° Ugh com P ar isons to existing sets of relevant 
data. Task II concludes with an identification of those aspects of thlwake 

analvtical nro^ requlr J ng furth * r development. Specific experimental and/or 
analytical programs to improve the model are recommended. 

The following sections describe and summarize the details of the rotor 

“un/tha^d^ "°f] deVel T' nt a " d reS “ US ° f the -t»di« ««l«d Z C 

published uudei' descrl f" :I ° n » f program and user’s manual Is 

pubn shed under a separate cover entitled "Development of a Rotor Wake/Vortex 

‘° e » ° ume User's Manual For Computer Program," (NASA CR-174850). 
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literature survey and data screening 




m 


1* General Electric in-house literature search syste 
2. MASA/RECON search system 
3* DoD/DTOC search system 

S^J^SSir 

that the data presented in the^“rLLnted 

‘ThT^ 1 ’' DetaUed — *lso availabU 

test JlV nr r 6 ; however » slnce it was intended to be one of the 
test cases for evaluating the developed rotor wake and vortex model these 

data were not used in the development of the empirical rotor wake Lei 
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Table I. Summary of Rotor Wake Reports. 
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Table I. Summary of Rotor Wake Reports (Concluded). 









3.0 EMPIRICAL PREDICTION MODEL FOR ROTOR WAKE AND VORTEX 


The measured stator upvash velocity spectrum exhibits tonal and broadband 
chacteristics . The tonal character of the spectrum can be ascribed to the 
fluctuating lift due to the periodic incidence of the rotor wake and vortex 
flow on the stator. The broadband character of the upwash velocity spectrum 
can be ascribed to the impingement of the turbulence of the rotor exit flow on 
the stator. The quasi-three-dimensional mod-1 of the rotor exit flow field, 
detailed in this section, employs a streamline-by-streamline approach. It pre- 
dicts the development of the rotor blade wake and secondary flow phenomena 
behind the rotor blade row. This in turn yields the stator incident gust and 
the corresponding tonal content of the stator upwash velocity spectrum, which 
has a direct bearing on the rotor/stator interaction noise level. An axisym- 
netric turbulence velocity spectrum model, described briefly in this section, 
predicts the broadband portion of the stator upwash velocity spectrum. 

This section briefly describes the methodology developed to predict the 
following: 


1. The viscous wake mean flow properties; namely, the streamwise vari- 
ation of the wake centerline velocity defect, the semi wake width, 
and the tangential wake profile 

2. The viscous wake turbulence properties; namely, the streamwise vari- 
ation of the peak turbulence velocities, and their tangential dis- 
tribution 

3. The secondary f low resulting from tip and hub vortices 

4. The stator upwash gust harmonic spectrum and the axi symmetric turb- 
ulent velocity spectrum 

Volume II of this report, "User’s Manual for Computer Program" (NASA CR- 
174850), gives a detailed description of the evolved computer program, a list- 
ing of the program, definitions of input/output parameters, and a sample 
input /output case . 


3.1 VISCOUS WAKE MODEL 


For the development of the empirical correlations of the wake center line 
velocity defect, semiwake width, and turbulent velocities, the extensive data 
reported in References 2 and 3 were employed. The semiwake width (6) is 
defined as the width of the wake at which the velocity defect equals one-half 
the wake centerline velocity defect (see Figure 1). 
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Isolated Airfoil Wake Rotor Wake 
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Figure 




3 * 1 * 1 Correlation of the Mean Flow Properties 


Empirical correlations employed in References 2 and 3 use different 
expressions for near-wake" and "far-wake" regions Regions with streamvise 
distance/rotor a.rodynaulc chord (s/c> lass than 0.5 have been designated as 

-z™ r :^: s : h ° se “ Ith s/c sreater than °- 5 *•" : 


An alternative form 
defect and semiwake width 
the course of this study, 
type: 


of correlation for both the wake centerline velocity 
using linear rational functions was developed during 
The correlations developed are of the following 


y - ax 1 b 

3 cx + 1 

where a, b, and c are empirically determined constants, and y could be 

W dc 5 $ 

^ or j or - and x = s/c 


( 1 ) 


ii??ntV ati( \ nal fu “ Ctl ° nS ° f thC type ShOWD in 1 are essentially 

inite series and are ideally suited for monotonically varying data such 

as wake centerline velocity defect and semiwake width, "ihe advantage Tftbe 

correlations using linear rational functions over the correlations developed 

the ent[re CCS \ “Vf® that onl y one expression is needed to correlate over 
the entire region of interest, and it is not necessary to provide a criterion 
to identify "near-wake" and "far-wake" regions. criterion 

The wake decay characteristics of five rotors [see Reynolds^) 

Si;::!' uwfii: haV e l„ studied. 

These data generally indicate that the higher loadings give larger wake veloc- 
ity defects. In References 2 and 3, local section d?ag coefficient Cn bis 
width U f 6d 48 3 cor f elatin 8 Parameter for wake velocity defect and semiSake 

tion ffr r a r 0 t r S K W ^ / ifferent loadin 8 s * As in the classical similarity solu- 
tion for a wake behind a two-dimensional isolated body( 10 ), the wake velocity 

to Tin \T Wi S th /°, r r ° t0r W3keS haVe been P-P° sad to be proportional 

defe!? ^ Z T °Z References 2 and 3 indicate that the wake^elocity 

ffse of f < T f ! Width f ° r r ° tors have a weaker dependence on C D than 
those of the isolated body. Also, in the case of a rotor, as a result of the 

semiwake°wi d t h ofs^^h th % adjacant blades * there is an upper limit to the 
semiwake width of S/2, where S is the blade spacing, irrespective of C D . 

hladff C ! rre !ci i 0 “ J OT the semlwake wi dth (5) has been developed with bo^ 
blade spacing (S) and rotor aerodynamic chord (c) as the normalizing dimension 

e or rotors with low solidity, rotor aerodynamic chord is probably the 


8 


C D in the ^-coordl nate* (such A , f r “ tlonal ° f of 

« were d -P e ?oT^a p r nts 

deviation, defined as (standard d / D the basis of normalized standard 

criterion for the seLcUon of rt evia ‘ ioa/ “ axi “- ^coordinate) x 100. The 
value of the normalized standard deviation? 1 C ° rrelation form was the minimum 

functions for ( s/s/verLs^/c)* (fi/s) f ^ a -2- btained with linear rational 
(s/ c)C d 1/8, respectively. ’ K * 1//C D versus (s/c), and 6/S versus 

The standard^eviation^of 6 the 

t'tiVT a sr 0 ; f 

Is the way the wake correlations « x P° nent * 8 1° the y-coordinate. This 

^ factor was ba«d on tha“a2 00 !^? ^ ‘a” 2 and 3. The 

(Reference 10). Now the standard j f f ° f tW ° dimens i° na l isolated bodies 

oaxlnun value of the y-coordlnate uTmH Th^ ‘“Of 1 " ls °- 530131 and the 

deviation of 13.361. I„ the case of £ le * d8 » uonnalltcd standard 

x-coordinate. The curvefit vIpIHc =. * l J* the C D exponent is in the 

maximum value of the v-ooord? 1 ^ standard deviation of 0.020405 and the 

a normalized standard deviation of^^Z 8 Whe 65 ^ Fisure 2) * Ttlls yields 
it can be noted that the correlation*!^ When comparing Figures 2, 3, and 4, 

better than that show^ L Figures 2 andT 4 fitS the data much 

functions for ( d/O^er^us^s/c)* (?/ V ? f ^ s -2- btained "1 th linear rational 
(s/O C D l/8, respectively? ^ <#/C> U ^ Versus (s/c) * and «/c versus 


fits^'datTo^rheite^tL” Lar:ho™ 4 ;„ t F i s :"r5 a "r6 8h °™ in Fis ? re 7 

atlon for (s/cTand ( «/S)? f 1 toe C ootes C ?h^t l> th e ° t | nor " alIzed standard devl- 
devfanon occurs If (,/sj and < Jc ) ar e^rr^ ^h' 

the linear rational function correlations of (a/^ i / Jr j C 9 * » 

with (a/c) vlaid largo errors. For t“a case of ' /s orV/”a (S/C) l/ ^D 
(fi/c) at large values of (c/n\ ^ ^ ° T correlated with 

(6/c) and («/S) will approach a constant ?al?I fun ^ tlon corr elation for 

thus Indicates that the semiwake width £ l 1 * in ^ e P endent of C D . This study 
equations developed for thf^miw^^dth a'e " deP6ndS WCakly ° n ^ 1116 
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Normalized Semiwake Width, 6/s 
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Figure 2. Correlation of 6/S with s/c for Midspan Data of Reynolds (Reference 
2) and Ravindranath (Reference 3) Rotors. 
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Normalized Semiwake Width, (6 /S) ( 1/ /Cj } ) 
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Figure 3. Correlation of (6/S)(l/v^) with s/c for Midspan Data of Reynolds 
(Reference 2) and Ravindranath (Reference 3) Rotors. 


Normalized Semiwake Width, 6 /S 



Figure 4. 


Correlation of 6/S with (s/c)(r 1 / 8^ c w., 

(Raf.rence 2) and Ravlndranath (Reference n “ ° f 


Normalized Semiwake Width, 6/e 



Figure 5. 


Correlation of 5/c with s/c for Midspan Data of Reynolds 
2) and Ravindranath (Reference 3) Rotors. 


( Reference 


lized Semiwake Width, (<5/c) (1//C|)) 



Figure 6. Correlation cf (o/c)(l/ /Cp) with s/c for Midspan Data of Reynolds 
(Reference 2) and Ravindranath (Reference 3) Rotors. 
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Normalized Semiwaie Width, <S/e 



Figure 7. 


Correlation of 5/c with (s/c)(C D i/8 ) for Midspan Data of Revnold 
(Reference 2) and Ravindranath (Reference 3) Rotors. 









(6/S) 


0.31875 (s/c) CqI/ 8 + 0.048 
0.268125 (s/c) C d ^/8 ♦ i. 0 


( 2 ) 


and 


u/ ^ . 0-2375 (s /c) - D l/S + 0.034125 
0-357 (s/c) c d 1/8 7 l.o 


(3) 


where 


5 = Semiwake width 

c = Rotor aerodynamic chord 
^ = Blade to blade spacing 

s = Streamwise distance from rotor trailing edge 
C D = Section drag coefficient 


« l SbL i ;.if i: x. e >??v.z:rr- te 

rotor blades. sxgmfies merging of wakes from adjacent 


Similar attempts to correlate the wake 
fractional powers of C D were made. Figures 
obtauied using linear rational functions for 
1//C D versus (s/c), and (W dc /W 0 ) (l/C D lA) 


centerline defect using different 
9, 10, and 11 show the curvefits 
( W dc/ w o) versus (s/c), (W dc /W 0 ) 
versus (s/c), respectively. 


11.28™™“::* t 11 * r i 

C"; (°? S , ta '' d " d showing that the If ^ 

Note however 0 thit^he (S/C) ^ lel< ? S the minimum normalized standard deviation. 
«ote, however, that the variation m normalized standard deviation for various 

CO e exp on ents i. -uch smaller than for the semiwake width correlation^ ?T g - 
is given bj: ° rreUtl ° n yieldln 8 the minimum normalized standard deviation 


( s 0-3675 (s/c) + 1.95 

V W o/\C D l/4 / /.bb (s/c) + 1.0 (4) 


where 

w dc = Wake centerline defect of the total relative velocity 

and 

W Q = Free-stream total relative velocity. 
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Total Relative Velocity Center 
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Figure 9. Correlation of W dc /W Q with s/c for Midspan Data of Reynolds 
(Reference 2) and Ravindranath (Reference 3) Rotors* 
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Figure 11. Correlation of ( W d c /W 0 ) (l/C^l/^) with s/c for Midspan Data of 
Reynolds (Reference 2) and Ravindranath (Reference 3) Rotors. 
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Thus, the velocity defect for rotor wakes is also seen to depend weakly 
on C D compared to the velocity defect of wakes of isolated bodies. 

The rotor leading manifests itself in terms of different free stream 
velocities and unequal wake widths on the pressure and suction sides of the 
rotor blade row (see Figure 1 for a conceptual illustration). References 2 
and 3 attempted to study the influence of rotor loading on rotor wakes in 
terms of the unequal wake widths on the pressure and suction sides of the 
rotor blade row. Figure 13 shows the variation of the rotor wake asymmetry 
(in terms of the difference in the wake widths on the pressure and suction 
sides normalized by the total semiwake width) and the semiwake width with 
streamwise distance for the data of Reference 3 at four rotor blade incidence 
angles. Similarly, Figure 14 shows the streamwise variation of the rotor wake 
asymmetry and the semiwake width for the data of Reference 2. Note that at 
some streamwise locations the suction side wake width is larger than the pres- 
sure side and at some other streamwise locations the pressure side wake width 
is larger. However, the extent of asymmetry in the wake widths decreases as 
the streamwise distance increases. This indicates that the difference in the 
free stream velocity on the pressure and suction side of the rotor blade row 
may be a more easily discernible effect of rotor loading on wake asymmetry 
than the distribution of unequal wake widths on either side of the rotor blade. 

The (S/S) dependence on (s/c) in Figures 13 and 14 seems to show different 
trends. In Figure 13, one is limited to (s/c) values of 3 0.6 wherein the 
near-wake effects are predominant, while in Figure 14, both near- and far-wake 

effects are evident. However, (5/S) values at comparable (s/c) locations are 
about the same . 

The wake model developed incorporates a constant inviscid velocity gradi- 
ent from the suction to the pressure side (see Figure 15). 4 constant inviscid 

velocity gradient in the absence of wake represents a linear variation of the 
velocity from the suction to the pressure side. This, when coupled with the 
wake model, simulates wake asymmetry on the pressure and suction sides, and the 
effect of rotor wake asymmetry on the stator gust description can be studied. 

In principle, the inviscid velocity gradient can be related to the rotor load- 
ing. As the rotor loading becomes greater, so does the difference between the 
free stream velocities on the suction and pressure sides, thereby increasing the 
inviscid velocity gradient. 

Unfortunately an empirical relationship between the rotor loading and the 
inviscid velocity gradient could not be developed in this program because 
applicable data were lacking. However, certain parametric studies were con- 
ducted (see Subsection 5.2.6) to enable a study of the influence of rotor wake 
asymmetry. This was accomplished by ascribing a specific value to the inviscid 
velocity gradient and studying the influence on the stator gust description. 

In the absence of any other specific information, the streamwise decay of the 
inviscid velocity gradient is modeled to be the same as that of the wake cen- 
terline defect. 

Next, using both a Gaussian distribution function and a hyperbolic secant 
function, the tangential distribution of the velocity defect of rotor wal >s 
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Normalized Streaowise Distance, s/c 

Figure 14. Rotor Wake Asymmetry and Semiwake 
Width Variation with Streatnwise 
Distance . 
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Velocity Velocity 




Suction 


Pressure Suction 


Pressure 


Rotor Rotation Direction 

a. Free Stream Velocity Variation Without Wake Defect 



Rotor Rotation Direction 
b. Free Stream Velocity Variation With Wake Defect 


Figure 15. Schematic of the Effect of a Constant Inviscid Velocity Gradient on 
the Tangential Distribution of the Free Stream Velocity With and 
Without Wake Defect. 




was correlated. To make a comparative evaluation of the Gaussian distribution 
function and the hyperbolic secant function, the tangential distribution of 
the wake defect of axial and tangential velocity components, reported in Ref- 
erences 2 and 3 at various operating conditions, was digitized and stored as 
a data base. 


Figure 16 shows the correlation of the tangential distribution of the 
normalized axial velocity component wake defect data with the Gaussian distri- 
bution function, and Figure 17 shows the correlation of the same data with the 
hyperbolic secant function. The Gaussian distribution function is given by: 

U d -ln2*n 2 

* e (5) 

dc 

where 

u d * Axial component of the wake defect at y tangential distance away 
from the wake centerline 

U£ C * Axial component of the centerline wake defect 
T) * Normalized tangential distance * x/(a/2) 

The hyperbolic secant function is given by 



sech (an) 


where a = cosh"* (2) * 1.3169579 


( 6 ) 


Figure 18 shows the correlation of the tangential distribution of the normal- 
ized tangential velocity component wake defect data with the Gaussian distri- 
bution function, and Figure 19 shows the correlation of the same data with the 
hyperbolic secant function. 

Note that the standard deviation of the curve fit for the Gaussian func- 
tion is slightly higher than that obtained with the hyperbolic secant function 
for the tangential velocity component of the wake defect (see Figures 18 and 
19), while the reverse is true for the axial velocity component. Therefore, 
both tangential wake profile correlations are equally applicable for describ- 
ing the tangential distribution of the wake defect. 


3. 1.2 Correlation of the Turbulent Flow Properties 

The motive for correlating the rotor wake turbulent velocities is to 
develop a turbulent velocity spectrum prediction model. The turbulence data 
reported in References 2 and 3 are in terms of the decay of the peak axial, 
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Normalized Axial Velocity Def 



Normalized Axial Velocity Defect, u^/u 
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Figure 17. Tangential Distribution of Normalized Axial Velocity Component Wake 
Defect Correlated with the Hyperbolic Secant Function 
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Normalized Tangential Velocity Defect, 



Figure 18. 


Tangential Distribution of Normalized Tangential 
Wake Defect Correlated with Gaussian Distribution 


Velocity Component 
Function . 
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Figure 19. 


Tangential Distribution of Normalized Ta - 
Wake Defect Correlated with the Hyperbolic 


ential Velocity Component 
- ecant Function. 


tangential, and radial turbulent intensities (turbulent velocity/local mean 
flow velocity). The peak turbulent intensities have been measured at the 
wake centerline where the mean flow velocity is the lowest. The turbulent 
velocity spectrum prediction model (to be discussed in Subsection 3.3.2) 
requires turbulent velocities and not turbulent intensities as input, so the 
turbulent velocity data of References 2 and 3 were renormalized with respect 
to a constant free stream velocity (velocity outside the wake region). For 
example, for yixial turbulent velocity, u'. 



where, 

u' * Axial turbulent velocity at wake centerline 

W c » Total relative velocity at wake centerline 

The first term on the right-hand side of Equation 7 is the axial turbulerr 
intensity reported in References 2 and 3. By using the correlation for the 
streamwise decay of the wake centerline velocity defect (see Equation 4) and 
the reported data on wake centerline turbulent intensity decay in the stream- 
wise direction (References 2 and 3}- the variation of the axial, tangential, 
and radial turbulent velocities normalized by the free stream total relative 
velocity at the wake centerline in the streamwise direction was obtained. 
Examination of the renormalized turbulent velocity data (see Figure 20) 
reveals that the axial and tangential turbulent velocity values are about the 
same, whereas the radial turbulent velocity values are higher. Thus, the axial 
and tangential turbulent velocities are correlated with one expression and the 
radial turbulent velocity is correlated with another expression. 

Linear rational function correlations were developed for the turbulent 
velocities, as in the case of wake centerline velocity defect data (see Sub- 
section 3.1.1), with different positive powers of Cp in the x-eoordinate [for 
example, (s/c)Cpl*^] and different negative powers of Cp in the y-coordinate 
[for example, (w' /W 0 ) (1/Cp^* 5) ] . These correlations were compared on the 
basis of normalized standard deviation, defined as (standard deviation/ maximum 
y-coordinate) x 100. Figure 21 shows the influence of the Cp exponent on the 
normalized standard deviation for the combined axial and tangential turbulent 
velocity correlation and Figure 22 shows the influence of the C D exponent on 
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standard deviation for tha radial turbulent velocity correla- 
. , , /u . 1 umu “ dbfnallted standard deviation occurs when (u’/W rt ), (v'/w > 

c^rilaM 05 *" “ tre }f ted “ lth (s/c)C D l-5. The linear ratlonal f2 nat i on 0> ’ 

(u'/U l or fv°/L re i f> ? i >8 t Y'“ norm alited standard deviation for 

/W 0 j or (v /W 0 ) is given by 


(“’/Vo) 
or / 

( V ’/W 0 ) 


425.0 * s /c * c 1 * 5 + 0.18 
, D 

12500 * s/c * C 1,5 + l.o 


( 8 ) 


i^ U r. 23 f Sh ° W ! th ! &0 ° d C ° lla P 8e of the axlal and tangential turbulent veloc- 

tnUtVl 5? r a ?° Ve COrrelatloQ * The linea r rational function correlation 
corresponding to the minimum normalized standard deviation for (w’ /W ) is 
given by: v 'o' 



o 


345.0 * s/c * Cp 1 * 5 + 0.264 

8625.0 * s/c * C D 1,5 + 1.0 


(9) 


Figure 24 shows the good collapse of the radial turbulent velocity data for 
the above correlation. 


The extent of anisotropy of the turbulent velocity scales in the rotor 
wake can be estimated by comparing the radial turbulent velocity variation in 
™L S fr direction to that of either axial or tangential turbulent 

velocity variation using Equations 8 and 9. Figure 25 shows the streamwise 
ariation of the ratio of radial turbulent velocity to either axial or tan- 
gential turbulent velocity. The anisotropy is high near the rotor trailing 
edge and decreases as the streamwise distance from the rotor trailing edge 
increases. However, even at large streamwise distance from the rotor trailing 
edge, a residual anisotropy exists. 8 


Hf . distribution of the normalized turbulent velocity is sig- 

nificantly different from that of the turbulent intensity, since the normalized 
turbulent velocity is obtained by dividing the turbulent velocity by the 

bv n dividino e rh 8 ^ ea J V * l0city > whereas the turbulent intensity is obtained 
by dividing the turbulent velocity by the local velocity. Reynolds (Refer- 

LTr iw T 1 “ dran f h (Reference 3) have shown that the axLl, tangential, 
K4uf dlal ^ tUrbulent lnteasit y components normalized suitably (see below) 
exhibit a Gaussian tangential distribution. For example, consider the tan- 
gential distribution of the axial turbulent intensity! Reynolds and 
Ravindranath showed that: 


T * T x.o -ln2n 2 

<-x> c ~ t x>0 G (10) 
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Normalized Radial Turbulent Velocity 
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and 

i 3 x/(«/2) 

x 3 Tangential distance from the wake centerline. 
By definition of the turbulent intensity, 


t « u ' (x) 
x U{7T 


(11a) 


T 

x,o 



and 


T 

x,c 


u* I x=o 


w r 


(lib) 


(11c) 


where 


u'(x) is the turbulent velocity at x distance from the wake centerline, 

W(x) is the local mean flow relative velocity at x distance from the wake 
centerline, 

w c ®ean flow relative velocity at the wake centerline. 

The data indicate that the free stream turbulence given by Equation lib is 
about 0.04 (see References 2 and 3). Substituting Equation 11 into Equation 10 
and solving for u'(x)/W(x), one obtains: 


u* (x) 
Wtx) 



0.04 


+ 0.04 


Cons ider. 


u' I x=0 

Wc 


u' I x=0 

Wo ' w dc 


u' lx=0 1 

W ° * ", 

~ w 0 _ 


( 12 ) 


(13) 
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where W dc is the velocity defect at the wake centerline. Also, consider the 
left-hand side of Equation 12, 


u' (x) _ u' (x) ^ w o 

W ( x ) W Q W( x ) 


Substituting Equations 13 and 14 into Equation 12 and solving for u'(x)/W n , 
yields 


\ W(x) 
u (x) _ 

Wo W Q L e 


-ln 2 h ; u'|x *0 


W„ 


0.04 | + 0.04 1 

f J 


(15) 


Let us look at 


W(x) 

Wo ' 


W(x) _ w o " w d(jO = x _ W d(x) = x _ w d(x) w dc 


W„ 


Wdc W D 


(16) 


where W d ( x ) is the wake defect at a tangential distance x from the wake center- 
line . 


The wake defect tangential distribution W d (x)/W dc has been shown to cor- 
relate well with the Gaussian distribution function (References 2 and 3), that 
is, 


W 2 

d(x) _ -ln 2 fT 

w dc 6 


Substituting Equations 16 and 17 into Equation 15, one gets 



We have developed linear rational function correlations for W dc /W Q (see Equa- 
tion 4) and — 7 - (see Equation 8 ). Denoting W dc /W 0 = f and u ® 


W„ 


W„ 


= g 


and simplifying, we obtain from Equation 18 
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u' (x) -ln2n 


W, 


- f e 


-ln4r| 


1-f 

2 r 


- 0.04 - 0.04 f 


- 0.04 


1-f 


+ 0.04 

The limiting forms of Equation 19 are: 
a. In the free stream, 

Lim u'(x) Lim 

- ■ M 

X > aj W/-v n -+<x> 


( 19 ) 


[right-hand side of Equation 19] * 0.04 


(20a) 


b. At the wake centerline. 


u*(x) 


x*0 


(20b) 


Thus, Equation 19 yields the expected values at the wake centerline and in 
the free stream. To study the tangential distribution of the turbulent veloc- 
ity, consider a numerical example. Let s/c - 0.5 and C D = 0.02. Equation 4 
yields 


«dc 


W 


o 


f - 0.1663 


Equation 8 yields 


u* |x»0 

W 0 


g a 0.0418 


As indicated by the limiting terms of the normalized turbulent velocity varia- 
tion in the tangential direction (see Equation 20), 


0.04 0.0418. 

w o 


Thus, we see that when the turbulent velocity is normalized by the free 
stream velocity the tangential distribution is quite flat, whereas when the 
turbulent velocity is normalized by the local velocity, it exhibits a Gaussian 
tangential distribution. The numerical example chosen is a typical case. For 
larger values of (s/c), the tangential profile of the turbulent velocity will 


41 


be even flatter. Hence, for the empirical model developed, which is intended 

the rotor wake characteristics for 0.5 to 5.0 chord spacings between 
the rotor and stator, it is reasonable to assume a flat tangential profile 
(that is, a homogeneous tangential distribution) for the axial, tangential, 
and radial turbulent velocities normalized by the free stream velocity. 


3.2 SECONDARY FLOW VORTEX MODEL 


The existence of secondary flow vortices in and behind a rotor blade pas- 
sage has been observed in experiments lo be dominated by the tip clearance 
leakage flow vortex formation, as discussed by Johnsen, Reference 11, and 
Lakshminaraj'ana , Reference 12. Other sources of secondary flow vortex forma- 
tions are present in rotor blade passages. These include end wall boundary 
layer secondary flow formations at both hub and tip, as discussed by Adkins 
and Smith in Reference 13. 

A tip clearance secondary flow vortex model was developed at General 
Electric for use in predicting secondary flow effects on fan tone noise, Ref- 
erence 14. This model was extended during the course of this study to include 
a hub vottex model to compute the velocity induced by both hub and tip vor- 
tices. The model consists of a forced vortex core and surrounding free vortex 
for both tip and hub vortirec 2S illustrated in Figure 26. A description of 
the prediction of the velocity fields created by the tip and hub vortices is 
given in the following paragraphs. 


3.2.1 Tip Vortex Model 

The tip vortex model superimposes a row of combined forced/free vortices 
and an image row opposite the annulus wall to enforce the zero radial velocity 
at the wall (see Figure 26). The tip vortices are located at a periodic tan- 
gential spacing equal to the blade spacing at the tip, S t . The tangential 
distance of the center of the vortex core from the projected downstream blade 
wake centerline is specified by the parameter, b t . The tangential location 
of the tip vortex has a significant influence on the stator upwash gust har- 
monic spectra, as will be shown through parametric studies in Section 5.0. 

From the experimental data surveyed, it is not possible to pin down a spe- 
cific value for b^, since it would vary not only from rotor to rotor, but 
also for the same rotor at different operating conditions. Hence, one will 
have to be aware of the sensitivity of the gust harmonic spectra to b t and 
excercise a judgment in the selection of b^. The blade span, tip clearance, 
end the tip vortex core size are designated by h, t, and a^, respectively. 

It is assumed that the vortex core contains all the shed vorticity so 
that the motion outside is irrotational . The radius of the tip vortex core, 
a t» determined from Rain's (Reference 15) analysis of tip leakage flow, 
which assumes inviscid flow and an idealized blade pressure distribution. 
Rain's expression to compute the radius of the vortex core at the trailing 
edge of the rotor blade is 







0.85 


f c 1/2 1 
-TE 0 ‘ 14 [t- <V tip J 


tipi (2D 

where T is the clearance gap height and (C L ) is the lift coeffici. 
the tip. C1 P 


. ent at 


rh c * rculatioti at the ci P Per unit span is related to (C L ) . through 

the following expresion: L tip 6 


(r) 


cip 


2 ‘Vtip - ( v tip * <c> ti P 


( 22 ) 


where (W 0 ) is the free stream velocity and (c) tip is the rotor aero- 
dynamic chord at the tip. 

The fraction of the tip circulation lost to the tip vortex (k)n was 
determined by Lakshmlnarayana in Reference 12 from blade pressure mealuraaents 
at the blade tip of a cascade rig. It can be approximated by 


(k) t ip = 0.23 + 7.45 (-|— ) for 0.01 < < 0.1 


(23) 


Thus, the amount of circulation in the tip vortex is given by: 

(rtip) vtx * OOtip (r)tip (24) 

The spanwise and the tangential <+y and +x directions respectively, see Figure 
6) velocities outside the tip vortex core were determined by employing an § 
extension of Lamb's (Reference 16) solution for the induced flow field about 
two infinite rows of vortices of finite radius as modeled in Figure 26. The 
resulting expressions for the induced velocities outside the domain R t follow. 


For the tangential velocity W t | 


tip : 


W tl 


(rtip)y tx f sinh M t 


tip 


T^L 


sinh N. 


cosh — cos p£ cosh - cos p^ 


(25) 


and, for the radial velocity W r | 


ti p : 


W 


(rtip) 


rl ti P 


2S 


vtx sin 

t cosh M t - 


sin p t 
N t - cos p t 


cos p t cosh 


( 26 ) 
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where 


M 


t 



T) 


N 


t 


2* 

— (y + a + T ) 
b t t 


and 



( 27 ) 

( 28 ) 


( 29 ) 


The parameters x and y are coordinates of the unwrapped annulus. 

The region within the domain R t behaves as a forced vortex (solid body 
of revolution). The domain is given by: 

(x-b t ) 2 + (y-a t -T)2 < a( .2 ( 30 


The angular velocity of the tip vortex can be related to the circulation of 
the tip vortex by: 


or = 


(Fti P ) vtx 
2* a t 2 


( 31 ) 


For a point P(x,y) within the domain R t , the tangential component of the 
velocity induced by the tip vortex is given by: 


w t ! 


tip 


- u 


r s m 


( 32 ) 


and the radial component by 



r cos ° 


( 33 ) 


where 


r = ((x-b t ) 2 + (y - a t 


t)2) 


1/2 
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and 


„ -1 f y - a t - T 1 

o * tan | — 

L x - b t J 


(34) 


; he t lndu H ed veio,:ity £ieid *» * -<*«i 
core. The radius of the vo«2 **! ?“* ° the presence of the tip vortex 

of the vortex in the streamvise directionT^d ^ ^ Vel ° City de ^iency 
dovmstr earn and the vorticity diffuses outwards Z^th 3 the ^ proceeds 
at any particular radius within th* ^ a ’ Also, the rotational velocity 

the velocity is negligible and the vortL “ reases until * ver V far downstream, 
mum strength of tX vort» £r. ?LT« .?*! C °“ Pletel S’ <*“>?«!• The maxl- 
to occur at the rotor blade tralllna .8 C i/** 8C * eaauise velocity i s assumed 
terminates. The d«a, . h ^ 8 “ h “* the flt >“ deaeration process 

viscous decay of a^tLd, u, ^ f' co " ™» after the 

of an airfoil Investigated bv Doean'ra' 18 .- 011 ? 1 trailln S vortex behind the tip 
velocity dlstrnu“onf l^th« -f J i R n efer '” oe 17 • »■ various 

■ ent with a linearized axisvmm.rrr. easured by Dosanjh were in good agree- 
nodel from Neuman (Reference 18) However^rh ncoopr esalble, viscous vortex 
length scales ' the 4K ‘ y ° £ the velocity and 

faster than predated b ? epproxlmately eight to ten times 

matic viscosity v . Since the vortex from S |-h hi ^ W * S based on a la minar kine- 
diffusion is stronger than for the nurelv i 6 . tip is turbul ent, the vorticity 
tual or apparent kfneJJl" c t L rnT “'““^ly, the vlr- 

greater than the l-mr ' lghb “«■ 

decay correlation, presented below, a constant 'v t / v Is assVed". V<>tt “ 

wise veloclty'wj^of "“f bUde CO “" tbe »«laum defect in stream- 

12) as: 7 dV the V ° rt “ COre ls Slveu by Ukshminar.yana (Reference 


w 


dv 


(u>a t ) ^ 
2U 1 


In 




0-75 cW m j 


(35) 


where 


i^ th rK main ! tr 630 fl ° W velocit y at the rotor inlet and 
'City through the hi ad* rv=.e.c,„„ and 


W m is the meai 


J. j. x yw veiocuv 

flow velocity through the blade passage. Also as th* u,c * 1 

stream, the maximum velocity deficit decays Is’ <1 V p ^ oceeds 

«: hS::r e Ti: 7 e-rL o i r „ro e fihe e o£ r d ~* y ° £ ■ 



(s/c) 




(s/c) 


) 


(36) 


Not e: 
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It is assumed that the tip 
edge. Hence, the decay or 
functions of (1 + s /c) and 


vortex originates at the rotor blade leading 
growth rates associated with tip vortex are 
not (s/c). 


Dosanjh et al. (Reference 17) have also shown the maximum rotational velocity 
due to the vortex decays with the downstream distance as (1 + s/c) (see 
note). It is assumed that the streamvise decay of the circulation of the vor 
tex follows that of the maximum rotational velocity. Also, at a given radius 
from the vortex core, the streamvise velocity deficiency is shown by Dosanjh 
et al. (Reference 17) to be: 


2 

W d (r , )/W < j v - exp [-In 2 (r'/a t ) ] 

where r' is the radial distance from the vortex core. 

From the above considerations, the three-dimensional velocity field 
due to the presence of the tip vortex can be calculated. 


3.2.2 Hub Vortex Model 

Though the existence of a hub vortex in rotors has been postulated, data 
presented in References 2, 3, 19, and 20 do not unequivocally confirm the 
presence of a hub vortex. Ravindranath (Reference 3) reports that the wake 
behavior near the hub region is considerably influenced by the secondary flow 
due to hub-wall boundary layer and also radial transport of mass and momentum. 

surveyed llteretur e «ef ereuces 2, 3, 5, 8, 9, 19, end 20 do ... to. eny 
kind of data which can be used to develop empirical expressions for hub vortex 
parameters such as vortex strength, vortex radius, etc.). Thus, only the 
methodology of calculating the velocities induced by a hub vortex is described 
in this subsection. When experimental data regarding hub vortex parameters 
become available, the emprirical relations can be firmly established and 
appropriate calculations of the induced velocities due to hub vortex can be 

performed. 

The mathematical development for computing the tangential and radial 
velocity components induced by the hub vortex follows a methodology identical 
to the tip vortex case. The origin and the axes of the coordinate system are 
still the same as for the tip vortex case. The changes relative to the tip 
vortex case are as follows (see Figure 26): 


Paramet er 

Tip Vortex 

Hub Vortex 

Radius of the Vortex 

a t 

a h 

y-Coordinate of the Image Plane 

0 

r ann “ r n 

y-Coordinate of the Vortex Center 

at + t 

r ann” r n~ a h 

Clearance Between Blade and Tip Wall 

T 1 

0 

x-Coordinate of the Vortex Center 

*t 

bh 

Tangential Distance Between Two 


s h 

Adjacent Vortices 

S t 

Circulation of the Vortex 

( r tip>vtx 

( r hub^ vtx 
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Since the identical coordinate system and sign conventions are applied 
tor both the hub and tip vortices, relations for computing the tangential and 
radial velocities induced by the hub vortex can be obtained from Equations 25 
through 34 (keeping in mind the differences between the hub and tip vortex 
systems listed in the previous table). Hence, the following expressions for 
the induced velocities outside the domain R h result. For the tangential 
velocity Wf- 1 , • 

w hub 


W 


(r hub ) 


vt x Slf * 

cosh Mft - c 


hub 2 

and, for che radial velocity W r 


sin 


os 


Ph cosh N h - cos ph 


] 


(38) 


W 


^ r hub 1 


vtx 


r hub 


2 S h 


[ 


hub ‘ 


sin 


Ph 


sin p. 


cosh M h - cos ph cosh N h - cos ph 


] 


(39) 


where, 

\ 




ann 


r, ) - (r 
h ann 


r h - v } 


(40) 


"h" y - <r a„n - r h> + <r a„„ * 'h * «„>} 


(41) 


and 


2x 

p h * <x - b h> 


(42) 


As in the case of the tip vortex, the domain Rh behaves as a forced 
vortex. The domain Rh is given by 


(x - b ) + (y - r + r. + a.) Z < a 2 

h ann h h - a h 


(43) 


The angular velocity of the hub vortex is related to the circulation of the 
hub vortex by 


(r hub) 


“h 


vtx 


2* a h ‘ 


(44) 
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For a point P(x,y) within the domain R h . The tangential component 
velocity induced by the hub vortex is given by 


of the 


W t! h ub “ ~ “h r sin a 


(45) 


and the radial component by 


W 


r ^ hub 


r cos a 


(46) 


where 


r 


[(x - b h ) 2 ♦ (y-r ann ♦ r h + a h ) 2 ] 1/2 


and 


tan 


-1 


/ ~ r ann * r h » a h 
x - b h 


(47) 


'•■'"rr 10 " 8 EOr ) Wl 1 *h a " d b h “»ld "on b e developed due 
to lack Of data on hub vortex, the folYSSing pseudorelations are incorpo- 
rated into the computer program with the goal of checking out the program 
logic for any compilation and execution errors! 


(r hub } vt x " °’ 2 * T (C L ) u..u * < W «>. . * (c) 


hub 


hub 


hub 


a h = °- 2 (c) hub 


(48) 


b h * 0*5 * S h 


where 


^ c ^hub * Rotor aerodynamic chord at hub 
^ C L^hub * Lift coefficient at hub 
^ w o^hub * Free stream velocity at hub. 


Once better definitions of the preceding parameters become available, the 
relations (Equation 48) can be replaced in the computer program. 


The model used for the streamwise growth of the radius of the 
and the viscous decay of the circulation, and the streamwise defect 


vortex core 
of the tip 
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vortex has been used for the hub vortex case, in the absence of anything 
better. 

The mathematical development for the hub and tip vortices neglects the 
influence of the annulus wall on the hub vortex and the influence o f the 
hub wall on the tip vortex. Also, since the wake/vcrtex model employs a 
streamline-by-streamline approach, only the tangential velocity component 
induced by the tip and hub vortices is computed. 


3.3 SPECTRAL PREDICTION OF THE STATOR UPWASH 


Sections 3.1 and 3.2 described the procedure to predict the mean and 
turbulent velocities of the rotor wake and the mean velocity induced by the 
hub and tip vortices in a frame of reference rotating with the rotor. This 
section briefly describes the following: 

a. Coordinate transformation of the combined flow field properties of 
the wake and vortex to a stationary reference frame, 

b. Evaluation of the description of upwash gust impinging on the stator 
and the corresponding evaluation of the gust harmonic spectrum, 

c. Evaluation of the turbulent velocity spectrum. 


3.3.1 Stator Upwash Gust Harmonic Amplitude Distribution 

The streamwise and the tangential velocity components of the rotor wake 
and the vortex flow in the frame of reference rotating with the rotor are 
summed algebraically to yield the integrated rotor wake/vortex streamwise and 
tangential velocity components. Thus, the wake/vortex integrated relative 
velocity vector is known. Adding the wheel veLocity vector to the wake/vortex 
integrated relative velocity vector yields the absolute velocity vector. Fig- 
ure 27 shows the velocity triangles in the free stream and at the wake center- 
line. For simplicity, it is assumed that the relative velocity has only a 
streamwise component. The perturbation velocity vector (V f ) in the station- 
ary frame of reference is obtained by the vector subtraction of the absolute 

velocity vector in the free stream (V 0 ) from that in the wake (V). The com- 
ponent of V* in the direction of V Q is the streamwise perturbation component 

and the component of V 1 normal to V Q is the upwash component. As the wake 
of one rotor blade passes in front of a stator, the upwash component of the 
perturbation velocity goes from a zero value (in the free stream) to a maximum 
(at the wake centerline) and to a zero value, again (in the free stream on the 
other side of the blade). From two-dimens ionaL nonstationary airfoil theory, 
the fluctuating lift has been proven to be acting at the 1/4 chord point. 
Therefore, the stator upwash velocities are evaluated at the 1/4 chord point 
of the stator. The fluctuating lift on the stator blade is one of the prin- 
cipal sources of the rotor/stator interaction noise. This subsection briefly 
describes the procedure employed for determining the spectral composition of 
the upwash component. 
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Notes 

_j* : Relative Velocity Vector at Wake Centerline 

Relative Velocity Vector in the Free Stream 
U: Wheel Speed Vector 

Ji: Absolute Velocity Vector at Wake Centerline 

V c : Absolute Velocity in the Free Stream 

V’: Perturbation Velocity Vector 

EB: Upwash Component of V ' v,J 

DE: Streamwise Component of V' V' 

y S 



Figure 27. 


Schematic Velocity Triangles for Rotor Exit Flow 
s 1 1 earn and at the Wake Centerline. 


in the Free- 
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The coordinate transformation yields the tangential distribution of the 
upwash component v£(y) over one blade spacing (S). Since V'(y) is periodic 
over S, it can be represented by a complex Fourier series as follows: 


vft(y) 



m 


imiry/(S/2) 

c m e 


(49) 


where 

i - /-f 

c m * ^ + ib m 

a^ b m are re *l nunbers. 


The above representation includes both the sine and cosine terms of the 

Fourier series. The complex coefficients c m can be evaluated using the ortho- 
gonality property of e «»*y/(S/2) function and ar<? given by 


S 


S 

f V’(y) e -i®»y/(S/2) 
o 


dy. 


• • • *“ 2 , — 1 , 0 , 1 , 2 , 


(50) 


The zero' th coefficient, c 0 , is essentially the average value of the upwash 
velocity over the blade spacing and as such of no interest in the spectral 
distribution, since it corresponds to the level at zero frequency. For a real 
function V n (y), * c m where c m is the complex conjugate of c m . Thus, in 
determining the gust harmonic amplitude, the coefficients of negative values 
of m do not contribute any new information. Therefore, the Fourier coeffi- 
cients corresponding to positive (nonzero) values of m only will be evaluated 
for this study. 

The integral in Equation 50 is evaluated nunerically in complex arith- 
metic by using the Simpson's rule for positive values of m (see Reference 20). 
The absolute value of c,„ and its phase are given, respectively, by 

1/2 

I c ml * (am 2 + b m 2 ) \ 


and I 

> m = 1,2, ... a. (51) 

' tan '‘ tv) 
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Note that | c m | has units of velo<-*ty. A nondimens ional decibel value for 
the coefficients of the Fourier series is defined by: 


dB * 10 log. - 
m *10 


'm' 


v ref 2 


, m - 1,2 . .. 


where V re £ is a chosen reference velocity. 


(52) 


While numerically integrating Equation 50 for the case of rotor wake only 
v without tip or hub vortex), it was noticed that dB,,, did not fall monotonical ly 
f ° r 1 “” ea81n 8 values of m, as has been experimentally observed, particularly 
when 4 /S approached a value of 0.5. The principal reason for this anomalous 
spectral behavior is that so far the influence of the wake from the adjacent 
blades has been neglected (see Figure 28). This is a reasonable assumption 
when «/S is small. However, when 4/s approaches 0.5, there is a large wake 
defect at the aidpassage. Figure 28 shows the case which has maximum defect 
at mid passage. The wake defect profile, whose Fourier components were deter- 
mined by the above process (that is, neglecting the wake from the adjacent 
blade), exhibits a sharp corner at y = 0.5 S (see Figure 28, wake defects for 
Blade 1 and Blade 2). A sharp corner in the profile excites higher harmonics 
(Fourier series for a step function). This was circumvented by combining the 
wa *® <? efects f rom the two adjacent blades. The combined wake defect does not 
exhibit sharp corner (Figure 28). The upwash spectrum of the combined wakes 
from the adjacent blades showed a monotonic fall with the harmonic number. 


The twist of the rotor blade from hub to tip causes a time/phase lag 
between the wakes from the hub and tip regions impinging on the stator. The 
analysis for predicting the phase lag due to the twist of the rotor blades 
was eveloped under a previous GE Independent Research and Development (IR&D) 
project and is included herein, for the sake of completeness. 


, . L ^ t the angle from the rotor blade stacking axis to the trailing edge be 
♦ Q (r) (see Figure 29a) and let the angle from the blade trailing edge to the 
stator leading edge along the wake centerline be A*(r,AX). It is assumed that 
the wake centerline "sheet" is fixed in space in a coordinate system fixed to 
the rotor. The angular displacement of the wake at the 1/4 chord point of the 
stator relative to the rotor stacking axis is given by: 


♦(r,X) - ♦o(r) + A*( r ,AX) 

The hub and tip values are respectively, 
♦h( r h» x h> = *o^ r h ) + ^(r h ,AX h ) 
♦ t (r t , x t) * ♦o^t) + & *(r t ,AX t ) 


(53) 


(54) 
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a. End View of Rotor Trailing Edge Line and 
Rotor Wake Line at Stator 1/4 Chord Point 


Rotor TE Stator 1/4 Chord Point 



Figure 29. Geometry for Aerodynamic Phase Lag Computation 
Due to Twist of Rotor Blade. 


Along a cylindrical stream surface, A* ls given by (see Figure 29): 


A«K(r, AX) = ~ tan 0 2 (r) 


( 55 ) 


where 0 2 <r) is the rotor relative air angle at the rotor trailing edge The 
^.railing edge angular displacement <fr Q is given by: 8 8 ' 


a / \ _ c(r) 

*o (r) ~ ~zT~ Sln y r 


(56) 


where c(r) is the rotor aerodynamic chord and y R is the rotor stagger angle. 

statorto^feel^h^ 5 he 1/4 chord P° int at a Siven radial location of the 

v k L Lt L ut "i * lmPBCt relative to the Phase at which the same rotor 
wake impacts the 1/4 chord point at the tip of the stator is then given by 

<Kr,X)l a g * 4<r,X) - <J> t (r t ,X t > 

The angular spacing for the rotor is given by 


(57) 


4 a — 

V B B 


(58) 


where B is the number of rotor blades, 
tip can be defined as 


A normalized phase lag relative to the 


♦(r,X) 


lag 


norm 


*(r,X) Ug 
♦ B 


(59) 


The normalized phase lag at a given radial location is 
to tip at that location relative to the rotor angular 


the phase 
spacing. 


lag relative 


3,3,2 Axisy°metric Turbule nt Velocity Spectrum Prediction 

hr na ^ e /° tal Stat !T 8USt Uf>WaSh s P ectrum ls treated as a superposition of a 
broadband spectrum determined by the rotor exit turbulence ch^acteristics 

rnmf s P ectrum determined by the rotor exit mean velocity cir- 

cumferentiai profiles. This section discusses the broadband component of the 
spectrum associated with the rotor exit flow turbulence. 

In Section 3.1.2, it was shown that, except in the immediate vicinity of 
the rotor trailing edge, the rms turbulent velocity is spatially uni f-> ^ 

the tangent ial drreetron when nonealrred by free-iLeL v^U="y. 1 Tgood 

Can be " odeUd as homogeneous turbulence. The mathematrcal 
°de I s for turbulence-rotor interaction given in References 21-23 were 
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selected as the starting point for the present turbulence-stator interaction 
model and associated stator gust upwash broadband spectrum prediction proced- 
ure. Only a few samples of turbulent velocity spectrum measurements were 
uncovered in the Literature Survey (Section 2.0); due to insufficient data, an 
empirical correlation approach was not possible. A theoretical model approach 
was therefore selected, and (as will be shown later) the model was verified 
using the few sample data cases available. 

Following Reference 21, a model of convected, homogeneous turbulence is 
assumed, incident upon a stationary blade row. The turbulence is assumed to 
be homogeneous, but not (necessarily) isotropic, with two of the three compo- 
nent directions having the same turbulence intensity and integral length scale. 
A sketch of the flow diagram, blade geometry, and nomenclature is shown in 
Figure 30, in the (x,y)-plane or cascade plane. The turbulence is assumed to 
convect with the mean (free-stream) flow, having axial velocity U and tangen- 
tial velocity Vp, and flow angle * tan“l(V^/U). 

Denoting axial, tangential, and radial directions by x, y, and z, respec- 
tively, and fluid— fixed coordinates by (x^, y 0 , z Q ) , we have that 


X Q = X + Ut; y Q = y - v yt ; £ o = 2 

Note that x is positive in the upstream direction, y is positive in the direc- 
tion of rotor rotation, and U is positive in the negative x direction. 

Assuming the turbulence to be statistically stationary in a reference 
frame moving with the fluid (x p , y Q , z Q - coordinates), the turbulent veloc- 
ities in the three component directions can be expressed as an integral summa- 
tion of shear waves (see Equation 1 of Reference 21) as follows: 

(u,v,w) = 

In terms of stationary coordinates (x,y,z), this becomes 

(u,v,w) =ydZ u> v,w^> exp(j[k*r + (Ukj - Vpky)t]l (61) 

In Equations 60 and 61, (u,v,w) are the (x 0 ,y 0 > 2 o) ~ components of the turbu- 
lent velocity, and dZ u>vw (k) represents the spectral density functions cor- 
responding to (u,v,w), given by 




dZ u , v , w (!$) exp (j k * r 0 ) 


(60) 


dZ u (k') dZ v (k") - 6(k'-k") * uv (k’ > dV dV (62) 

The component of velocity normal to the stator chord line (see Figure 30) 
at an angle Y s relative to the axial (x) direction is given by 
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v n * u s in Y s + v cos y s 


( 63 ) 


Thus 


I v nl 


2 . 2 

u sin Y s + 2 uv s ^ n Y s cos Ys + V 2 cos : 


Ys 


(64) 


s r n”“* t r r r nt v \ iocity f °*- «ioc ity „ 

me, that is, the upwash or gust velocity spectrum Yn _, i 


.2 

♦nn ~ *uu Sln Y s + 2* uv sinY s cosy s + cos 2 y s 


normal to the 
Ynn> then given by 

(63) 


at zero ^incidence ^"^rftagMr'b^hl’* ^ * flat 1 pl f te stator < and rotor) 
stator and extends the turbulence model tre3tS 3 Staggered 

higher in complexity than th* • » ) , axis ynmetnc case, one step 

21. Equations 61 and 65 become identiCa/to^h^ot 0 * 161 aa8X f med f n Ref erence 
21 if Vy is replaced by -V - -*r thlttt J h « r t0r equations m Reference 
“Qj., the rotor relative air ancle a ° r T e6 ^ 8peed > and Ys is replaced by 

=55= sa= = : zL? 

appropriate equation., since References 22 .!/.,( v y and n r by - Ts ln the 
the rotor problem in Reference 2l"=L«£.:\T 1 l< ’ P * d e «' msi ‘>" s 
results and transfer. them dir eetly to the stator p“u«?“’ “ °* e th ° s ' 

blade C row^ e fhe turbulence S spectru. t for ^th tUr ^ U t*" C t i ? cident u P° n a rotor 
to the relative free-streL v^c“ y vectt *°™ 1 


. 2 


♦nn = ♦uu sin a r - 2 sin+ r cosa r * uv + ^ cos 2 ar 

For isotropic turbulence, the component spectra are given by 

♦ 1 _ 

4,k2 ) 6 ij k 2 j 

where 


(66) 


(67) 


* 2 ‘ k l * k 2 2 * k 3 


( 68 ) 
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and 


E(k) 


8 k 4 u 2 
r L (L“ 2 + k 2 )' 


( 69 ) 


Upon combining Equations 66 and 69 by first evaluating Equation 67 for ♦ = 

*11’ = *12 and *vv = ♦ 22 > then substituting the results into Equation^ 

whe following expression for is obtained: 


, E(k ) 
nn 4ir k 4 


k) f, 2 2 

T— k cos a + 2k k, sina cosa 
•+ l 1 r 12 r r 

2 A , 2 1 
n a + k.. 
r Jj 


+ k 2 sn 


(70) 


Following the format of Reference 21, we need to integrate Equation 70 over 
spanwise wavenumber k 3 * kj (note that axial wavenumber k= * ki and tangen- 
tial wavenumber kp * k 2 notation is adopted herein). Integrating * nn (kt k-i 
k 3 ) over -•<k 3 <- gives the gust spectrum ♦£ n (k 1 ,k 2 ) as follows: 


OD 

/ 


♦ nn ( k l» k 2^ * / ^nn^l »^2 * k 3^ dk 3 


(71) 


Substituting Equations 69 and 70 into the above, we have 


nn 


OD 

f 2 ^2 

(k coso + k_ sina 2 ) + 
x r y r z 

/ ir 2 L 

(L-2 + k_ 2 + k_ 2 + k 2 )3 
x y z 


dk_ 

z 


(72) 


In the preceding expressions, L denotes the turbulence integral length scale, 
and u 2 denotes the mean-square turbulence velocity. 

By employing the substitutions k 2 = A tan 8, where A 2 = L 2 + k- 2 + fc— 2 
Equation 72 can be integrated analytically to yield the foil — - * y ’ 


owing: 


- u 


r nn 


t 


4xLA' 5 L 3(k ^ C ° S “r + k y Sln<1 


r > 2 ** 2 ] 


(73) 
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After rearranging terras, this can be rewritten as follows: 


* 

♦ 

nn 


2 

u 

4irL 


j k x * c 2 k ; k ; * 


2 r-2 

+ L 


(L -2 + k_ 2 + k_2) 

X y 


5/2 


( 74 ) 


where 


C 1 1 + 3 cos^a r 


c 2 * 6 sma r cos a r 
c 3 * 1 + 3 sm^a r 


(75) 


The next step is to integrate 


k y • 


nn iver circumferential or tangential 
The turbulent velocity spectrum is a function of 
, S ~y* As 1S discussed in References 21 and 22 k~ (avi^i 

18 reLted to the spectru " b ^ rough tz: 


wavenumber k£ * 
both and k' 


w » 2irn V/D - !Jk x 


(76) 


IiatId n wi"J 1 discrete U ? ber C ° rrt \ Sp ° nds . £o the spanning raode lobe number asso- 
cludes -hit -h P i l frequency harmonic disturbances. Hence Reference 21 con- 
-«rnr . integration over k*. for the purpose of estimating turbulence- 

rotor interaction noise, need only be over the ran*e nf v- “ * Cbrbulence 

propagating or cut-on modes m the fan duct * y Corres P° bdin S 

Although the actual turbulence spectrum is the resiil*- of ► i * 

tion ri f U -h° n - ° Ver Sl \ Wavenumbers k y* ^ is of interest to evaluated por-^ 
tion of -he tangential uavenumber spectrum that car produce acoustic enerev 
Therefore, an rntegration of ♦£„ „ as carrred out over finite corre- 

>i“‘* »* Ui, duct at the frequenci^ 

Mach numbers of interest. The integration limits can then be extended 1^- 
to determine the total turbulence spectral amplitude at frequency i " 

evaluate 1 " 1 " 8 ** th * £urbulence 8 ust velocity spectrum, we need to 


* (k_) = 7^ 
n x 4irL 


*(k_) 

x 


(77) 
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where 


•GO 

x 



C I k x * C 2 k H k y * C 3 k ? 2 L ' 2 

( l “ 2 ♦ k 2 + k 2 ) 57 " 2 

* y 


dk_ 

y 


( 78 ) 


and a and b are the lower and upper 
dimensional) cutoff. For a rotor 
21 and 22) 


limits of ky corresponding 
these limits are given by ( 


to duct 
see Refe 


( t wo- 
ences 


a - - X/A-H 2 a - 2ir n /S 
b * + x//l-M 2 - 2w n /S 

3 


(79) 


where 


* Axial flow Mach number U/c Q 
X • Frequency wavenumber w/ c Q 
s - Rotor blade-to-blade spacing (pitch) 
n * An integer, < n < <» 


f tHat 8eference 21 had an erroneous exponent of 1/2 

To do *Ltl en y Equatlon 78 - Equation 78 can be integrate 
To do .his we again introduce the transformation- 


in the denominator 
analytical ly . 


A 


2 


1 + k_ 

x 


2 


and 


k_ 

y 


- A tan 0 


where 

ke * Lk and k_ * Lk_ 
x x y y 


Then 

dky = A sec^0 d0 


and Equation 78 can be written 


as 


2 r h 

* (k x } = Ju J j l 1 + k_ 2 ] cos 3 0 + [c 2 k. A] cos 2 e si 


n0 


0. 


+ (C3A 2 ] sin 2 0 cos0 


d6 
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This can be integrated analytically to give the 

♦(kj): 


following expression for 


"V - £ { [‘ ♦ ‘1 s * 1 ] [' 

~ (2 + cos 2 e ) sm 8 1 
a a 


(2 ♦ cos 0 ) sin8, 
b b 


where 


- [<= 2 ^-aJ j^cos 3 e b - cos 3 0 a J 
+ Jc 3 a2 J ^sm 3 0 b - sin 3 0 a J j 


and 


e a * tan" 1 (I/A) and 0 b - tan _1 (b/A) 


a = La - -(L/S) (xS//l-M + 2im) 

a 


b * Lb = + (L/S) (xS//Tm 2 - 2m) 


( 80 ) 

(81) 


(82) 


Equations 77 and 80 through 82 define the turbulent velocity upwash spectrum 
for isotropic turbulence incident upon a rotor. The fraction of upwash energy 
that produces propagating modes is accounted for by the limits I and b. 

Having outlined the methodology to derive the gust upwash turbulent veloc- 
ity spectrum using the simple case of isotropic turbulence incident upon a 
rotor, we now consider the case of axisymmetric turbulence incident upon a 
rotor. From Reference 22, the component three-dimensional wavenumber spectra 
analogous to Equations 67 through 69 are given by the following equations: 


(kj + 

k^) F(k lt k 2 ,k 3 ) 

(83) 

" k l k 2 

F(k 1 ,k 2 ,k 3 ) 

(84) 

(k* + 

^3 ^ ^(^1 >^ 2 *^ 3 ) +> k^ G(k^ ,k^) 

(85) 
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where F and G are given by 


2 

*2 


i r 4 2 

L L a 

at a 


(1 • lV * lV) 

a a t t 


2 L L ? 


2 u 2 - u 2 (1 * L 
t a 


x J 


<1 ♦ L 2 k 2 ♦ L 2 k 2 ) 
a a t t 


and 


2 2 
k +• lc 

2 K 3 


> k . 


In the above expressions, 

~2 

u a axial component of turbulence intensity 

T 

transverse component of turbulence intensity 
L a * axial turblence length scale 
L t = transverse turbulence length scale 


( 86 ) 


(87) 


( 88 ) 


It is assumed that the axis of symmetry of the turbulence is in the axial 
direction, and the turbulence intensities and length scales m the two compo- 
nent directions normai to the axis of symmetry are equal, being denoted by ^7 
c ’ r ® s P ec ^ lvel y- This allows a specification of radial turbulence mten- 
scll e and length SCale dlfferent from axial turbulence intensity and length 

EouatL T S “ bStltUte E ^ a tions 83 - 88 into Equation 66 and the result ir:to 
Equation 71. Carrying out the algebra and subsequent integral evaluations as 
was done for the isotropic turbulence case, we get the following expression for 
f nn- 


nn 


L L u 2 
a t a 

4* A 5 


tAJ * E ♦ e^kjL.k, * c,L 2 k 2 


(89) 
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where : 


1 * L a k l * L t k 2 


c i = E + 3 ujl ) 2 cos 2 a 

a r 

c 2 * 6 s ina cosa 

w 3 r r 

c 3 = E + 3 sin^a 

r 


( 90 ) 


sm 2 a r + e cos 2 a r 

‘ 2 - ‘ VV * 


For isotropic turbulence, u? - u ? and L .r 

Thus expression 90 reduces to „ r L; i L 9’ SO that 6,1 and E » 1 . 

ulence. Note thit L t /L a appears in cLh 6 ” Ved fonnul ae for isotropic turb- 

expressions for cj, c 2 , and^-,. " C ® b J n J tlon cos a r m the above 

in the two parameters E and (L*/l ) coT tt" 6006 axisymnetr y manifests itself 

c ' sl ' wos a r . 

expressions for^t^t^'turbu ^ f ° U ° Wln * 

b 


♦ Ck_) 

X 


/ 


E * C l k l * c 2*1*7 + c 3 k 2 ‘ 


5/2 


dk? 


(91) 


° + ^1 + * Z } 


formulation 7sJ excep? thatch* tera^^the h3S ^ **** f °™ ? S £he isot «-opic 
I.' ♦ C l*l 2 - The integral evaLatioThas III l .* c l^l 2 i« replaced by 

ication. We therefore have the following clc^ld-^ Wlth the . same term modi- 
symmetnc turbulence incident upwash gust^pLtrl oTa^r 8 : 10 " ^ ^ aX1 ’ 


*(kj) 


a 


30 ♦ k, 2 ) 


[ E * c l k l 2 ] [' 


(2 + cos 0, ) sin0. 

b b 


~ (2 + cos e ) sme 
a a 


[• 


3 3 

cos “ cos 9 
D a 


]‘[ C 2 k l <l * k i 2 )’ /2 ] 

] * [ C 3 * k V>] [ 


sin 6, - si 
b 


A] 


(92) 
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where 


^5 * tan _1 ( b//l + k x 2 


®a * tan 1 ( a /^1 + i ^ 2 j (93) 

and ' ' 

a ® L t a j b a L^b • 

It was previously shown that the stator inflow turbulence spectrum 
can be derived from the rotor inflow turbulence spectrum by replacing a r by 
-Y s . Equations 92 and 93 can thus be applied to stator inflow turbulence with 
c l> c 2 > an< * c 3 replaced by the following: 


c - E ♦ 3 (L /L ) 2 cos 2 Y 

1 t a s 

c- * -6 (L /L ) sinY cosY 

2 t a s s 

. ? 

C 3 ■ E + 3 sin Y s 

E - sin 2 Y + c cos 2 Y 
s s 


(94) 


The stator axisymmetric upwash turbulence spectrum is then given by: 

~2 

W ■ -TO" ,(k i> (,5 > 

a 


where *(ki) is given by Equations 92, 93, and 94. The limits a and b given by 
Equation 79 or 82 are simpler for a stator because the integer n results from 
a summation over blade-passing frequency harmonic numbers, and only the n «0 
term applies for homogeneous turbulence incident upon a stator blade row. Thus, 
a * -b, and this allows simplification of Equation 92 to the following form: 


2 L 


•(k^ 


_ 2 

3(1 + k L ) 


|[ Etc i\][ 


(2 + cos * ) sin®, 
b b 


■] 


+ 




2 

03d + ) 



(96) 


Equations 93 through 96 completely desc rib e the stator incident turbu- 
lence upwash spectrum, and require as input u u£, L a , and L t . This formu- 
lation is used in the computer program. 
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4.0 ANALYTICAL WAKE MODEL STUDIES 


Two analytical wake model studies were conducted during the course of 
this research program to complement the empirical wake model study described 
in Section 3.1. They were: 

!• Analytical Correlation of Turbulent Wake Decay by Dr. C. Hah of 
General Electric. 

2. Rotor Wake Momentum Integral Analysis by Prof. B. Lakshminarayana 
of the Pennsylvania State University. 

A brief description of these studies is given in the following sections. 

4.1 ANALYTICAL CORRELATION OF TURBULENT WAKE DECAY 

The objective of this study was to provide a physically sound formulation 
for the development of a semiempirical model for the decay of mean velocity 
and turbulence intensity downstream of a turbofan rotor. As the primary appli- 
cation of the program is for predicting noise and unsteady vane forces due to 
rotor— stator interaction, the analysis is focused on the moderately near and 
far wake region (0.5 to 5.0 chord lengths downstream of the trailing edge of 
the blade). 


4.1.1 Mean Velocity Decay Formulation 

The evolution of turbulent wake behind a turbofan is governed by various 
length and velocity scales. Depending on the dominant scales, the entire wake 
field can be divided into several regions (very near, near, and far wake). 
Because of these distinctly different scales, one universal correlation for 
the entire wake field may not be feasible. In the following, a formulation 
is developed th.it is valid in the moderately near and far wake regions (approx- 
imately 0.5 to 5 chord lengths downstream of the trailing edge.) 

The Far Wake - Far downstream of the trailing edge, the effects of load- 
ing, rotation, curvature, etc. on the wake evolution disappear and a similarity 
solution based on slow evolution can be derived. In this region, the differ- 
ence between free stream velocity (W Q ) and local velocity (W) is smaller than 
free stream velocity. 

w d - W Q - W << W Q (97) 

With the above assumption, the governing momentum conservation equation in the 
streamvise direction can be simplified to the following (Reference 10): 

a W d l 3T t 

~ W o 'Ss“ = P TS~ (98) 
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where 0 is density and t c is turbulent shear stress. In the tar wake region 
simple mixing length theory can be used and the decrease ot W d is assumed to’ 
be a power function of 3 . Also, the mixing length is assumed to be propor- 
tonal to the wake width, which is the spacing of the cascade (S) in the far 

wake. Then the wake defect (W d ) should decrease in proportion to s~l (See 
Reference 10)* v 

Schl ichting , using the above assumptions (Reference 10), obtained the fol- 
lowing analytical distribution of velocity behind a row of cylinders: 


u = W ° / s ¥ S / 2*n 

d i^U/s cos (— 


(99) 


where L is mixing length and the experimentally obtained value of L is 


s' « 0.103 for - * 0.125 (c » Chord) 

th ^ tur * ) °f an wake. Equation (99) should be approximately valid in the far 
wake where the effects of blade loading, machine rotation, and other operating 
parameters disappear. Note, however, that Equation (99) gives only a single 
harmonic amplitude for the spectrum of the wake waveform, which does not agree 
with the experimental evidence. 

_ • r f ? ea f “ ak g .~ To extend the above analysis to the moderately near wake 

region, the following general form of mean velocity decay was assumed: 


Me 


W„ 


= A + 


s+D 


( 100 ) 


where W dc is the maximum wake defect and D represents different scale effects 

n6ar ^ P ° SSlbU shlfc of origin of the wake. A can be used as 

another parameter for the wake shift (such as Mach number effect). However 
A is dropped for the present derivation. * 

effecr h nf f n rm I! la i t a O a (100) Sh ° Uld be correct at far downstream, where the 
• rhP c . should disappear. In the far downstream, the dominant parameter 
is the solidity of the cascade (c/S), and the following form is proposed for 


B 


- *. (!) 


( 101 ) 


°i l0 f? ing ’ raachine rotation, blade geometry, etc., should be 
included in D. The following simple form is assumed for the current study: 

D = A 2 C D 1 s P (102) 

where C D is blade drag coefficient and p>l . 
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Combining Equations (100), (101), and (102) yields: 


^dc 

Wo 


A! (S/c) 3 
"*! C D? S ’ P 


( 103 ) 


The formulation (103) can be reduced to the same form as the exact solution for 
the far wake given by Equation (99), and the coefficients A 2 and p for the near 
wake effect, along with A^, can be optimized with data. It should be noted, 
however, that Equation (103) yields W^ c * 0 in the limit as s * 0, which is 
physically unrealistic. Thus, selecting A^ and A 2 to give a satisfactory data 
match may prove difficult. 


4.1.2 Turbulence Intensity Decay Formulation 

For the general correlation of the turbulence quantities, a far wake solu* 
tion is first derived. Then the formulation is generalized for the moderately 
near wake region. 

To evaluate proper length scales for the decay of turbulence quantities 
at far wake, the following two scales are defined: 

L s * the sc^le of change in the streamline direction 

Ljj * the seal* of change in the normal direction. 

then 



and 

- uv * o(u 2 ), u 2 * o(y 2 ), v 2 = o(m 2 ) 

where U is turbulence velocity scale, u and v are fluctuating velocity compo 
nents, and 0( ) indicates the order of magnitude. 


From the continuity equation, an order-of-magnitude analysis yields the 
following: 


V - 0 



(105) 
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Using the order analysis of (104) and (105), the cross-stream momentum equation 
in the far wake can be simplified as follows (see Reference 24): 


3v^ 

3n 


P 3n 


(106) 


Integration of Equation (106) yields 
P . ~2 p o 

— + V * — 
p p 


(107) 


The streamwise momentum equation in the far wake can be written, 

„ 3U „ 3U 1 3P 3 T 3 Zv~ ^ vf3 2 U ^ 3 2 u\ 

3s + 3n P 3s 3s U 3n ^3 S 2 j n 2 J 


Combining (107) and (108) 


„ 3U „ 3v 3.2 2. 

U -5— + V -s— + -j— (u - v ) + 
3s en d s 




(108) 


(109) 


The orders of the magnitude of the terms in Equation (109) are estimated 
with Equations (104) and (105). For the high Reynolds number and in the far 
wake , 

— ■> - and » 0 
v L s 

Then the first and the fourth terms remain large and should balance each other. 
The order of magnitude of these terms is: 


„ 3U „2 w dc / w o U L n\ u 2 
3s " U L s “ ^ W * U ‘ L s ) L n 


( 110 ) 


and 


3n 

Tl re fore 


(uv) * - r- 


( 111 ) 


W. 


o ^ ^dc 
V U 


IT * ° (L n> 


( 112 ) 
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I W 


In the far wake , 



Therefore , 



0 (L ) 
n 


(113) 


(114) 


The derived Equation (114) indicate 
tne same order as the velocity defect, 
should also be proportional to the s“i, 


s that turbulence intensities are of 
The decay of turbulence intensities 
as the mean velocity decay. 





Ti * (S/c)3 
1 “T 3 

s + T 2 C D 2 s 


(115) 


where T 


l> t 2, T 3 , 


should be determined from the experimental data, 


becomes •* anisotro P^, but 

-r— 2 a 13 

be explicitly included in T 2 . § ’ ^ nUm ’ rotatlon » etc -) can also 

The correlation method used in inn \ i .* u 
relation function: * S ^ ase< * on c ^ e following cor- 


V = «LUl 
cx + d 

The above equation can be rewritten as follows: 


(116) 


y = A’ + 


b 1 

x+c ' 


(117) 

oirsZri:Zi" e °z “n r t ‘”«L d ;:? very " eu - uo " ev "- 

decay properly far downafreao ip i “ ies do pot 

correlations, however the formnlpH , S - * 0ver the range of data 

, u wcver, cne tormulation works quite well. 

dowpsSe::^: dased onthe .« «« 

effect of blade spacing is explicitly included, although 
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the dependency is somewhat arbitrary. 3y optimizing the coefficients of tne 
corrections these formulations can be successfully used both in the moder- 
ately near and far wake regions, any additional effects on the decav of rh.=> 

lif l J 1 C Chi y wioni U Seh.i?o^ “In f and 

Srr-S “ <" ~r 

(lib). It was therefore decided to use Equation (116). q 


4.2 


ROTOR WAKE MOMENTUM INTEGRAL ANALYSIS 


technia u r al r iS ° f 3 ^ W3S carried out using the momentum integral 

chnique. Two momentum integral equations on each side of the wake orovided 
the necessary closure equations to predict wake deferr f ^ 


^•2.1 In t roduc t ion 

turbulent wi th ^ubstant ial eSS ?^ ° r fan rocor blade is three-dimensional and 

depend on nations blade and'flov parameters t” d ^ 

blade loading, Reynolds number, «Lh numb"/ uje“ r lot ^ 

the earlier analyses available for the rotor wake are based ok “«1, “ 

orofT?e.° nS " — , a " d Laksh ” i "»rayana (Reference 25) assumed tnat velocit? 

5I?ec is"L n“ " “a U " aa ‘ i2ad tha equations (assuming that the ve ocity 
defect is small compared to the freestream velocity) to derive an.lvticiH 

i° : ° f the - ”- ai 


W = | “8(A+$)v^ + U 2 (BScosA) 2 ) 

Ue U e | ~ 4Bflsin\(s 2 cos 2 A)v T f 


where A and B are given by 

A = B 2 R s 2 and (A $)/B = ir 2 /4 



A^cos 


i n-*] i/2 " 


where u c and U e are strearawise velocity defect at 
stream velocity (streamwise) , respectively. W is 
and ♦ is a rotation parameter given by 


the centerline and free- 
the local radial velocity 


fl 2 S 2 

Ue 2 


2 n 

{sin Acos^A} 
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All other notations are explained in the nomenclature (Section 8.0). 

and t^r^r (Reference 26 > ^ed the momentum integral approach 

^nd IV equations to predict the decay of the axial, tangential, 

•rends in /? ° C1Cy ^°“P onent s, and they showed good agreement with measured 
flow Thi a * " Speed fan rotor ‘ ^ analysis is restricted to free vortex 

to other ^ a " d 

tn i R rr; y> K ~ l and HirSh (Reference 27 > used the Pennsylvania State data 
to indude the effect of streamwise pressure gradients in the decay of the 

an tL V ! lr = 0rrelation a 8 reed with the Pennsylvania State data 

on the rotor wake characteristics. 

Hah and Lakshminarayana (Reference 28) developed a computer program to 

stress modei aCt T^is a is° nS °f . m ° tion f in Chree dimensions, employing a Reynolds 
stress model. This is an elliptic formulation with initial conditions (at 

e*ori»I a ^ lln f 6dg ® derived from the experimental data. They also provided an 
expression for the near-wake characteristics for two-dimensional wa Tes 

The present analysis is a new formulation and the analysis is valid fr»r 

Tthlff “ rj- 38 C l! tar Wake * 1116 defect » velocitiel c^uld L large 
this formulation. The analysis and equations are presented in this report* 
results from the analysis will be evaluated at a later date P ’ 


4.2.2 


Momentum Integral Analysis 


is t m ° St Suitable coordinate system for the analysis of the rotor wake 
with the’/; r T °, m in FlgUre 31 ‘ The ^reamwise coordinate aligned 

Ts the r H /I! SC Y eamllne is ^ven by s, n is the principal norma/ an/r 
s the radial direction. The corresponding velocities are u, v, and w respec- 
tively. The analysis given below is valid for an asymmetric and incompressibl f 
rotor wake, even through its extension to include the compressibility effect 
is fairly straightforward. The profile in the n direction Is assume! toll 
Gaussian m the momentum integral formulation. 

given^Ln Re^ence coordinate system shoun in Figure 31 are 


Radial Momentum Equation : 

.. 3 W 3W 3u 

W 3r + V 3n + U 17 + 2 ° VcosX + 2«UsinX 

- - lip*. i_ (T j 

P 3r 3n ^ rn' 


u 2 . 2 V 2 2, 
— sin X - — — cos X 


where 

* r 2 

P = P o~ 


(118) 
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sure f'**' 

approxmat ion) , the foliouing approxhaation can be aide’ * Uy " 


i i£l „ 1 3p „ . Uf 2 

0 3r p 3r 2QU e siiU + — sin X. 

Hence, Equation (118) reduces to 

W I? + v f£ + U H - 2fl(U e - U)sinX 4 ^L-~ 1,2 
03 r 


2 a 

sin X * — ( T ) 
3n VT rn ; 


where r rn is the turbulent shear stress in the radial direction. 


(119) 


Streamwise Momentum Equation : 

If the streamwise pressure gradient (3p/3s) 
momentum equation reduces to (Reference 29). 


is neglected, 


the streamwise 


U -f^ ♦ V — 

3s 3n 


an 

+ W — - 2flWs inX + 


UW 

r 


. 2 , 

sm X 


3t 


sn 


3n 


( 120 ) 


Continuity Equation : 


3U 3V 3W W 

3s 3n + 3r r 


( 121 ) 

on the suction side, and a similar in f0 o„r' enter to the edge of the wake 
center to the edge on the pressure side* l °” ** carrled 0,,t llcm th « u * ke 

boundary ‘«««» “• ««"tu. integral analysis for a 

erence Jl>. LM Ind Ir Unrated t«o-di.en,ional flows (Ref- 

s/Jr - 0. The continuity Iguanirthln’lldlllrto 111 "'' 10 "’’ " ltt ' “ * 0 and 


3V 

3n 


3U 

3s 


( 122 ) 


iU 

T 3U 

/ 17 


dn + C . 


(123) 
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The velocity V is not zero at the wake centerline, so 


V = 



V 

c 


(124) 


The appearance of V c on the right side is the essential new element in tne 
momentum integral analysis for a wake. Using the standard procedure (Refer- 
ence 32), the following momentum integral equation car be derived for two- 
dimensional viscous flows: 

<* 8 1 _ a H ♦ 2 dU e _ V c / U e - U c ^ 
ds * 1 U e dx * U e l U e f 

where ® i is the monentum thickness, H is the shape factor, and U c Is velocity 
at the centerline. 


4.2.3 Momentum Integral Analysis for the Rotor Wake 

Following a procedure similar to the two-dimensional flow described 
above, the momentum integral equations for 3D flows can be derived by inte- 
grating Equations (119) and (120) in the normal (i») direction. The velocity 
V in these equations is derived from the continuity Equation (121), 


V 



(125) 


and substituted in Equations (119) and (120) before integrating. For the sake 
of brevity the detailed derivation is not given here, but is similar to the 
3D boundary layer momentum integral derivation given by Mager (Reference 33). 
The additional terms appearing in the wake equations arise because of nonzero 
values of V c in Equation (125). The momentum integral equations governing 
the wake on each side of the wake center is given by: 


30 11 . 2 3 U e a . 39 L2 . 2 3 U e a .1 3 U e , ! 


3 s Ue 3 s U 
sin 2 X 


0 . _ + — 


, 3n 

6 + ± ® 6 


U e 3r 12 U e 3s 1 U e 3 r 2 


(0 +0 ) - 2 
r ^ 12 IV 


^2^sinX 


/ Ue ~ Uc \ 
U e V U e / 


(126) 


and 


39 


21 + 2 _ 


3U 0 30 

0 __ + 


e a . ~”22 . 2 3 ^e a . sin 2 X 


3s U 3s 


21 


+ - 9__ + 


3r 


U 3r 


22 


(6 - 0 - 6 ) 
K 22 11 1 


+ 26 


ftsinX 

1 a. 


v c / ^e Wc 
= ^ \~^~e 


(127) 
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where 



( 128 ) 


(129) 


(130) 


(131) 



(132) 


where 6 is the distance beyond the wake edge 


and there are nine unknoJSs^e™* g 126) 6 3nd (127) for eac h side], 

side (total of fifteen knnwne »» t/ 12 ’^ 2l y ??’ f** 6 2» v c» u c > w c on each 


the wake). Hence, some simplifying^,-^ ’ ti„ V Sam “ f ° r boch sides of 
couple various thicknesses in the ^ C ° made in order to 

in selecting the velocity profiles U(n ‘and'uf ? f aSSUm P tion can be made 
;n low .peed teat, 1„ Rei.reeces 2 1 jl"'’: U has «U-.Mabl iahed 

both the St Teamwise and the rad a "foci 5 «e -»iia, for 

by velocity. The similarity profiles are given 


" 1 “ u c exp (-0.693 n 2 ) 

W w c 

* u~ exp ( _0 -693 n ) 


(133) 

(134) 


where 


n 



> 


or 


n 



S. L 


*s 


surfaces of ^he^^^respect ivel^ Pr6SSUre 3nd Suction 



Substituting Equation (133) and (134) 
can be expressed in terms of velocity 


in (128) through (132), the thicknesses 
defects. For example, 


'll 


■ L f !■ * u c exp(-0.t>93 n 2 )] lu c exp(-0.693 n 2 )l dn 
L ^[l.06 - 0.7530 


since 


flD 


-a 2 x2 /* 

dx ' IT 


Similarly, 


W r u~ 

• l 2 - 0.7530 


21 


22 


, r w c u c w c i 

= " l L»; 1 - 06 ■ ufuf 0 - 7530 J 


LW C 2 


[0.7530] 


u c 

L jf- 1.06 
u e 


"L u; 1*06 


(135) 


(136) 

(137) 

(138) 

(139) 

(140) 


Substituting Equations (135) through (140) in Equations (126) and (127), we 


s momentum: 


3F 2 / 3u e\ 3 / W c \ 2 3 U e / W \ 3 

3s + U e' 3 s ) f+ 3?(°- 71G u7) + u^T7-( °- 71 G u^) + !£ 

H 3G e ^ sin 2 * ( n 
~ U e 3r r (°* 

.. v c ( u e - u c ^ 

u 0 \ UZ ) 


3s 


u e \ w r 

71 G - F— / 
u c / U 


Of. \ W r 2Hflsin> 


( 141 ) 
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r momentum: 


3 

7? 


where 


(- F £) (' F £) ♦ i? (-0” “ £) * f e ^(-0.71 B £) 

( 142 

* ^ (-0.71 H £ - F - <A * 2 G 2i*“i . Ic / «e - w e \ 

\ U e ) V e U e \ U e j 


u e [ l - 


06 - 0.7530 


G - 1.06 L ~ 
u e 


(143) 


(144) 


H * 1.06 L 

°e 


(145) 


E « ch side of the wake has two equations: Equations (141) and (142) for 

a 1 nr r i vo unlrtArjn^. t t . .. ' * 


a totai of five unknowns; Lp, L s , u c , W c , and V c . The velocity U e ( s ,n)'is 
either prescribed or derived from an inviscid analysis. 

An additional closure equation is needed to complete the set of e^ua- 

I tor oredirtlna f" ho V9ri ef i ~ C- XT IT * . _ ^ 


tions for predicting the variation of V„ 

a. i j a _ C > 


‘ , . 7"~ — ’c> u c» w c» i*p » and L s with distance. 

related U fh^ nd A “f° lx (Refer « nce 31 > employ Head's entrainment equation. Head 
elated the growth of a turbulent boundary layer to the amount of fluid the 
boundary layer entrains from the irrotational flow around it (see Reference 
34), to yield an auxiliary equation which came to be known as the Head's 
entrainment equation. In the case of a wake at sufficient downstream dis- 
tances, the entrainment was shown to be proportional to the velocity defect 
at the wake center line (Reference 34). By definition, the mass flow rate 
within a boundary layer of thickness « is directly proportional to («-«,), 

th ” J fi C 6 lsp } a ^® ment thickness. Hence, the streamwise variation of 
the mass flow rate within a boundary layer (in other words, entrainment) is 

devIlon^d € ^h by rJ he 8 ^ eamwi - se Privative of («-«!). Coustieux and Aupoix have 
«hich U ,i«n bjr dl " e,,Sl0nal » E the Heed's entrainment equation 


i a r 1 36 . 

u ( 6-6 ) - : 

Ue 3 * L 6 1 J 3 r 


C E + — 

E U a 


( 146) 


where Cg * Head's entrainment factor. 
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4.2*4 Calculation Procedure 


On the pressure side, there are two equations governing the rotor wake 
flow (Equations (141) and (142) with Lp in place of L) . There are also two 
equations on the suction side (Equations (141) and (142) with L s in place of 
L) . These equations, along with Head's entrainment equation (146), can be 
solved simultaneously to derive the values of the wake width (L s + L D ). 
velocity defect (u c ), radial velocity at the wake center (W c ), and the wake 
trajectory (V c ) as a function of the downstream distance. These equations can 
be programmed to predict the wake properties. During the initial calculation, 
the terms 3U e /3s can be neglected. It may be possible to correct for the 
change in U e through a blockage coefficient derived from the wake width and 
profile. The programming and calculations were not carried out during the 
present study due to a lack of time and resources. The analysis does, however, 
offer promise of the ability to predict rotor wake structure, accounting for 
many more influential parameters with much less empiricism than the current 
correlation approach (Section 3.1). It is strongly recommended that this 
effort be pursued as a follow-on study to the current program* 
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5.0 MODEL EVALUATION 


The empirical rotor wake/vortex model developed in Section 3.0 was evalu- 
ated through data-theory comparisons for three rotors (Rotor 55, JT15D fan 
rotor, and Rotor 67) and various parametric studies. Table II lists certain 
key geometric and aerodynamic parameters for Rotor 55, JT15D fan rotor, Rotor 
67, and various parametric studies. Section 5.1 deals with the data-theory 
comparisons, and Section 5.2 deals with the parametric studies. 


5.1 DATA-TdEORY COMPARISONS 


The data-theory comparisons are performed for rotors which were not part 
of the data base used to develop the empirical wake relationships, so that this 
set of data-theory comparisons constitutes an independent validation of the 
empirical model. The inviscid velocity gradient, WTIV (see Figure 15), has 
been set equal to zero for all the data-theory comparisons. Subsection 5.1.1 
describes the data theory comparisons for Rotor 55; Subsection 5.1.2 addresses 
the same for JT15D fan rotor; and Subsection 5.1.3 deals with Rotor 67. 


5.1.1 Rotor 55 Data-Theory Comparisons 

The rotor wake and stator upwash harmonic data for Rotor 55 were obtained 
from References 19 and 6 for three rotor speeds. The geometric data and aero- 
dynamic data at the three rotor speeds were obtained from Reference 35. The 
section drag coefficients at the required span locations were computed from 
the measured profile loss coefficients, reported in Reference 35. 

Figure 32 compares the predicted and measured axial variation of center- 
line wake defect (W^ c /W 0 ) and semiwake width (6/S) at 30% span from tip for 
80% and 115% design rpm's. The data at 80% design rpm were taken with and 
without the tunnel velocity. The data at 115% design rpm were taken with a 
tunnel velocity of 41.1 m/s. The wake defect predictions at 80% design rpm 
agree reasonably well with the data. The theory overpredicts the seraiwake 
width for the 80% design rpm case. Both the predictions and data indicate that 
the influence of the tip speed at constant loading (that is, for operation 
along the fan operating line) on the wake characteristics is negligible. The 
agreement between the data and predictions for the semiwake width at 115% 
design rpm is very reasonable. 

Figure 33 shows the comparison of the measured and predicted gust harmonic 
amplitude and broadband turbulence spectra at 30% span from tip for the 80% 
design rpm case. The levels of the data are adjusted so that the prediction 
and the data for the BPF level of the gust harmonic amplitude agree. The mea- 
sured broadband level is adjusted by the same AdB and compared to the predic- 
tion. The measured and predicted gust harmonic spectra show similar trends 
(monotonically falling with increasing frequency) for this spanwise location, 
which is close to midspan location where effects of the tip vortex are not 
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Table II. Key Geometric and Aerodynamic Parameters for Rotor 55 
JT15D Fan, and Rotor 67. 


Parameter 

1 

Rotor-55 

JT15D Fan 
Rotor 

j Rotor-67^^ 

Diameter, Inches 

20.00 

21.00 

-4 — 

19.54 

Number of Rotor 
Blades 

15 

28 

< 

22 

Number of Stator 
Vanes 

25*2) 

66 

34 

Hub-Tip Ratio 

0.460 

0.434 

0.477 

Tip Clearance, Inch 

0.0315 

0.0236 

0.0157 

Design rpm 

8020 

15800 

16043 

Tip Speed at Design 
rpm, fps 

700 

1448 

! 

i 1407 

Rotor Solidity at Tip 

0.896 

1.34 

1.327 

Rotor Solidity at Midspan 

0.998 

1.53 

1.772 

Rotor Solidity at Hub 

1.186 

— 

2.879 

Camber at Tip, ° 

22.65 

N. A. * 3 ) 

11.41 

| Camber at Midspan, * 

39.54 

N.A.* 3 ) 

17.57 

Camber at Hub , * 

45.22 

N.A.* 3 ) 

54.12 

Aerodynamic Chord at Tip, Inches 

3.651 

3.086 

3.655 

Aerodynamic Chord at Pitch, Inches 

3.033 

2.536 

3.646 

Aerodynamic Chord at Hub , Inches 

2.412 

— 

3.720 

Drag Coefficient, Cn 
at Tip* 4 ) 

0.107 

0.075 

0.034 

C D at Pitch* 4 ) 

0.038 

0.031 

0.012 

C D at Hub* 4 ) 

0.128 

— 

0.102 
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Table II. Key Geometric and Aerodynamic Parameters for Rotor 55, 
JT15D Fan Rotor, and Rotor 67 (Concluded). 


Parameter 

Rotor 55 

JT15D Fan 
Rotor 

Rotor 67^^ 

Work Coefficient, We 
at Tip(*> 

0.486 

0.314 

0.382 

W 0 at Pitch^ 

0.688 

0.371 

0.509 

W e at Hub ^ 2 3 4 ^ 

1.098 

— 

1.095 

Rotor Outlet Relative 
Flow Angle, 82 at Tip * ^ 4 ' 

40.2 

63.1 

60.6 

&2 at Pitch, *( 4 ) 

19.1 

52.4 

41.8 

02 at Hub, *( 4 ) 

-4.9 

— 

- 6.7 


Notes 


1. Rotor 67 is a two-stage rotor. The parameters tabulated here 
correspond to the first stage. 

2. A 25-vane stator was employed in the wake study of Reference 19 and 
an 11 -vane stator was employed in the aerodynamic performance study 
of Reference 35. 

3. N.A. - Not available. 

A. Quoted Aerodynamic Performance for Rotor 35 is at 80% design rpm; 

Quoted Aerodynamic Performance for JT15D fan rotor is at 66 % design 
rpm; 

Quoted Aerodynamic Performance for Rotor 67 is at 100% design rpm. 
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Figure 33. Comparison of Measured and Predicted Gus c Harmonic Amplitude 
and Broadband Turbulence Spectra. 
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significant. The predicted gust harmonic spectrum shows steeper fall with fre- 
quency than the measurement, since the predicted semiwake width is larger than 
the measured semiwake width. As semiwake width increases, the wake profile 
over a blade-to-b lade spacing resembles a sinusoidal function with a period 
equal to the blade-to-blade spacing. Hence, the contribution of higher har- 
monics to the spectrum decreases sharply and the BPF dominates the spectrum. 

The band width for the broadband data and predictions is 50 Hz. The predicted 
broadband spectrum shows a steeper fall with frequency above 1.5 kHz compared 
to the data. However, they do show comparable trends. Keep in mind though, 
that the measured broadband data show a significant amount of scatter, and the 
best estimate of the random data was made for this comparison. 

Figures 34 and 35 compare the predicted and measured gust harmonic ampli- 
tude and broadband turbulence spectra at two spanwise locations (10% span from 
tip and 50% span from tip) for the 96% design rpm case. The agreement between 
the measured and predicted gust harmonic amplitude levels and the broadband 
spectra for the 10% span location from the tip is good. This spanwise location 
is influenced by the tip vortex. The tangential location of the tip vortex 
for this computation was selected to be b t /S - 0.42. The agreement between 
the data and predictions for the broadband turbulence spectra is also good. 

For the case of 50% span location from the tip (see Figure 35), the data show 
that the even multiples of BPF show higher levels than the odd multiples (for 
example, 2xBPF level is higher than BPF, 4xBPF level is higher than 3xBPF) . 

This is a typical trend if there are two defects within one blade-to-blade 
spacing, as would be the case when a strong tip vortex exists in that region. 
However, at this spanwise location, which is sufficiently far from the tip 
region, one would anticipate negligible influence of the tip vortex. The pre- 
dictions indicate a monotonic fall-off with increase in harmonic number, since 
the prediction model is based on the rationale described above. At present 
this behavior of the measured gust harmonic spectra is not understood. The 
predicted and measured broadband levels show reasonable agreement. 

Figures 36 and 37 compare the predicted and measured gust harmonic ampli- 
tude spectra at three axial spacings for the 80% design rpm case at two span- 
wise locations (10% and 30% span from tip). The tangential location of the 
tip vortex is chosen to be b t /S = 0.42, as before. For the case of 10% span 
from tip, both the data and theory predict the 2xBPF levels to be higher than 
BPF levels at all three axial spacings (Figure 36). As described earlier, 
such a trend occurs whenever a tip vortex is dominant. Both the data and 

theory indicate that the influence of the tip vortex is considerable, even at 

the farthest axial spacing (X/c * 1.77 at the pitch streamline). The trends 

predicted by the theory are quite representative of the data. Figure 37 com- 

pares the predictions and measurement of the gust harmonic spectra for the 80% 
design rpm case at a typical ^idspan location (30% span from tip). Note the 
monotonic fall-off of the harmonic content at the midspan location, showing 
that influence of tip vortex at the raidspan location at all axial spacings is 
negligible. At smaller axial spacings, the semiwake width is smaller. The 
wake shape resembles a delta function when the semiwake width goes to zero. 

In such a case, the Fourier coefficients of higher harmonics do *ot fall off 
rapidly. This situation is observed at X/c - 0.54. As the wake spreads out 
at larger spacings, the fall-off rate of the harmonics increases. The trends 


86 



UNCLAS 











Velocity 
Wheel Speed 


• Rotor 55 

• 96% Design RPM 

• 10% Span from Tip 

• X/c = 0.54 at Pitch Streamline 

• Inviscid Velocity Gradient, WTIV = 0 

• Band Width for Broad Band spectrum = 50 Hz 


5XBPF 


4XBPF 


Predicted Broad 
Band Spectrum 


Measured Broad 

Band Spectrum (Reference 19) 


Frequency, r'Hz 


Figure 34. Comparison of Measured and Predicted Gust Harmonic Amplitud 
and Broadband Turbulence Spectra. 




Rotor 55 
96% Design RPM 
50% Span from Tip 
X/c = 0. 54 at Pi tch 

Inviscid Velocity Gradient, TIV = 0 
Band Width for Broadband Sp^trum = 50Hz . 

7 ■ — - 


2XBPF 


+ XBPF 


o 


3XBPF 


Measured Gust 
Harmonic Amplitude 
(Reference 19) 




Predicted Gust Harmonic 
Amplitude 



Predicted Broadband Spectrum 


Measured 

Broadband Spectrum 
(Reference 19) 





Harmonic Number. 












predicted agree well with the data. The nu.:.imum deviation in data and predic- 
tion is seen for X/c = 1.23. The measured semiwake width at this location is 
smaller than the predicted seraiwake width (see Figure 32). Thus, the predic- 
tions at X/c s 1.23 show a steeper fall-off rate for higher harmonics than do 
the data. 

Figures 38 and 39 are similar comparisons of the data and predictions of 
the gust harmonic amplitude for the 115% design rpra case at two span locations. 
The data at 10% span from the tip indicate that a strong tip vortex (Figure 38) 
is absent. Hence, the predictions for this case do not include any tip vortex. 
Note that the agreement between the theory and predictions is good at X/c - 0.54 
and 1.23. The predictions at X/c - 1.77 show a similar fall-off rate as the 
data, except for the 4xBPF tone. Figure 39 compares the data * r d predictions 
at 30% span from the tip location. The agreement is quite good at X/c * 0.54 
and 1.7, and the fall-off rate of data and predictions agree with each other 
at X/c = 1.23. 

In summary, the empirical rotor wake and vortex model is seen to predict 
very well the observed trends both in the gust harmonic speeds and the 
turbulence spectra for Rotor 55. 


5.1.2 JT15D Fan Rotor Data-Theory Comparisons 

The relevant data for JT15D fan rotor data-theory comparisons were taken 
from Reference 36. The wake upwash component was measured at four spanwise 
stations on the bypass stator and the gust harmonic data at four speeds (43%, 
54%, 61%. and 66% design rpros) were supplied for the data-theory comparisons 
(Reference 36). Figure 40 compares the spanwise distribution of certain key 
geometric and aerodynamic parameters relevant to rotor wake studies for Rotor 
55 and for the JT15D fan rotor. The aerodynamic data for Rotor 55 are for the 
80% design rp® case, and the aerodynamic data for JT15D fan rotor are for the 
66% design rptn case. The section drag coefficients for Rotor 55 were computed 
from the measured profile loss coefficients. Vne section drag coefficients 
for the JT15D fan rotor were computed from the profile loss coefficient, which 
in turn has been correlated with the diffusion factor for NACA 65-( A^Q)-series 
and double-circular-arc blades (Reference II). The axial spacing for Rotor 55 
corresponds to the case when X/c at the pitch streamline equals 0.54 (Note the 
high rotor solidity and the large rotor exit relative flow angle of the JT15D 
fan rotor compared to those of Rotor 55). 

The rotor wake and vortex computations were performed with the tip vortex 
tangential location of b t /S * 0.45. The gust harmonic amplitudes are computed 
using a 10 fps velocity as a reference velocity so that a one-to-one compari- 
son to the data could be made. Figure 41 compares the measured and predicted 
gust harmonic amplitude spectra at four spanwise locations of the bypass sta- 
tor at 66% design rpm. The predicted gust harmonic amplitude spectra at all 
the spanwise locations are in very good agreement with the data- In order to 
evaluate the spanwise extent of influence of the tip vortex on the gust harmonic 
spectra, computations were performed by suppressing the tip vortex calcula- 
tions. Figure 42 compares the measured and predicted gust harmonic spectra 
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Figure 39. 






























with and without tip vortex computations at the four spanwise locations for 
the 66% design rpm case. The axial spacing between the rotor and the bypass 
stator of JT15D fan, normalized by the rotor aerodynamic chord, is approxi- 
mately 0.70 (Reference 36). For the above rotor-stator spacing, we note from 
tj-gire 42 that the tip vortex influences only the 11% span location from the 
tip and has no influence at the other spanwise locations. 

Figures 43, 44, and 45 compare the .data and predictions of the gust har- 
monic amplitude spectra for 61%, 54%, and 43% design rpra's respectively. For 
all these calculations, the tangential location of tip vortex is the same as 
that for the 66% design rpm case, namely, b t /S - 0.45. The data and predic- 
tions agree very well at all the spanwise locations at 61% design rpm (Figure 
43). The agreement between data and predictions for the 54% design rpm case 
at 11%, 57%, and 90% span of the bypass stator is excellent (Figure 44). Data 
and predictions at the 33% span locations disagree on absolute levels, but 
show similar fall-off rates. For the case of 43% design rpm, the data at 
the tip streamline (11% span) show a more pronounced tip vortex contribution 
than is predicted by the model- The agreement can be improved by adjusting 
the tangential location of the tip vortex. But, owing to the excellent 
agreement between data and predictions at three out of four speeds, the tip 
vortex location was maintained at b t /S - 0.45 for this speed also. Hie data 
and predictions agree well at the 57% and 90% span locations, and show similar 
fall-off rates at the 33% span location. 

In summary, the agreement between the predicted and measured gust har- 
monic amplitude spectra for the JT15D fan rotor at various span locations for 
various rpms along the engine operating line has been found to be excellent. 


5.1.3 Rotor 67 Data-Theory Comparisons 

Rotor wake measurements made on Rotor 67 at NASA Lewis using a Laser 
Doppler Velocimeter were supplied to us (Reference 36). This formed the third 
data set for evaluating the empirical wake model. The rotor wake measurements 
supplied were at the peak efficiency point, and essentially included blade 
passage averaged relative velocity wake profile data at various spanwise and 
axial stations. Rotor 67 is a two-stage fan having low-aspect-ratio first 
stage blading. The rotor wake measurements supplied were for the first stage 
rotor. These rotor wake measurements were obtained with the first stage 
rotor-alone configuration. The geometric and aerodynamic data for Rotor 67 
were obtained from Reference 37. As in the case of Rotor 55, the section drag 
coefficients were computed from the measured profile loss coefficients. 

Wake centerline defect normalized by the free stream velocity and semi- 
wake width normalized by the blade-to-blade spacing data were extracted from 
the wake profiles of the relative velocity data at the peak efficiency point 
(100% design rpm). At this operating point, a considerable portion of the 
first stage rotor has supersonic wheel speeds. Measurements revealed that for 
about 50% of span from the tip, the shock loss coefficient is equal to or 
greater than the profile loss coefficient, implying that more than half the 
span cf the rotor blade has supersonic flow over it (Reference 37). Only the 
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profile loss coefficient was considered in computing the section drag coeffi- 
cient. Using the above information, rotor wake predictions were made. 

Figure 46 compares the measured and predicted wake centerline defect and 
semiwake width axial variation at 50% span from the tip. The agreement between 
tie predicted and measured semiwake width axial variation is very good. How- 
ever, the measured centerline wake defect does not agree with the predictions 
for X/c < 0.4 Rotor 67 has a significant amount of supersonic flow over the 
blades, and as such may have strong shock-boundary layer interaction which 
could complicate the flow character at the rotor blade trailing edge. This in 
turn could result in a wake decay vhich is quite different from the wake decay 
of the low speed rotors used in developing the prediction model. Thus, the 
measured wake centerline defect does not agree with the predictions at dis- 
tances close to the trailing edge. However, for X/c > 0.4, the measurements 
and predictions agree well, presumably because the complex flow character 
present at the trailing edge has dissipated and the wake displays the typical 
characteristics of wakes of subsonic tip speed rotors. 

Figure 47 shows a comparison of the measured and predicted wake center- 
line defect and semiwake width axial variation at 10% span from the tip. As 
in the case of 50% span from the tip, the predicted wake centerline defect 
agrees with the measured wake centerline defect for X/c > 0.4. However, the 
semiwake width predict ons do not agree well with the data. At 10% span from 
the tip location, the tip vortex significantly affects the rotor wake shape, 
and as such the measured wake at this spanwise location is highly complex in 
shape. Hence, extracting semiwake width data from such complex wake shapes 
involves a considerable amount of uncertainty. The experimental data of semi- 
wake width at this spanwise location show almost no variation in the semiwake 
width over the range of axial stations. This is ascribed to the tip vortex 
interactions with the rotor wake. 

In summary , the predicted and measured wake centerline defect in the mid- 
span location for Stage 1 of Rotor 67 are in reasonable agreement for X/c > 

0.4. The agreement for the semiwake width in the m: dspan location is good. 

The predictions and data of the wake centerline defe t do not agree close to 
the rotor trailing edge. 


5.2 PARAMETRIC STUDIES 


The reasonably good agreement between the theoretical predictions and 
the experimental data shown in Section 5.1 for three rotors gives credence to 
the developed rotor wake and vortex model. The next step in the evaluation of 
the model is to perform parametric studies to determine the influence of vari- 
ous geometric and aerodynamic parameters on the rotor wake and vortex charac- 
teristics. The parametric studies included the following: 

1. Influence of rotor-stator spacing on the rotor wake characteristics 
and the stator upwash gust harmonic spectrum. 
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Figure 46. Comparison of the Measured and Predicted Wake Centerline 
Defect and Semiwake Width Variation at Midspan. 
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2. Spanwi se variation of rotor characteristics and their influence on 
the wake characteristics and the stator upwash gust harmonic spec- 
trum. 

3. Influence of rotor solidity on the rotor wake char ac t er is t ic s and 
the stator upwash gust harmonic spectrum. 

4. Influence of section drag coefficient on the rotor wake character- 
istics and the stator upwash gust harmonic spectrum. 

3. Influence of rotor loading at constant tip speed on the rotor wake 
characteristics and the stator upwash gust harmonic spectrum. 

6. Influence of the inviscid velocity gradient on the rotor wake char- 
acteristics and the stator upwash gust harmonic spectrum. 

7. Influence of the tip vortex parameters (tip clearance and tangent al 
location of the tip vortex) on the stater upwash gust harmonic spe * 
trum. 

8. Influence of the axial and transverse velocity and length scales on 
the broadband turbulent velocity spectrum. 

5.2.1 Influence of Rotor-Stator Spacing 

Rotor-stator spacing is an important parameter in the rotor-stator inter- 
action studies. As the rotor-stator spacing increases, the rotor wake dissi- 
pates, reducing the fluctuating lift force component on the stator and thus 
lowering rotor-stator interaction noise. 

Figure 48 shows the influence of axial spacing on the normalized wake 
centerline defect and the semiwake width. The computations performed refer to 
the 10% streamline of Rotor 55 at 80% design rpm. The pertinent aerodynamic 
and geometric parameters were taken from Reference 35. The wake centerline 
defect is about .6% of the free stream velocity at X/c =0.5 and decreases to 
about 4.5% at X/c = 5.0. The normalized semiwake width (w/S) reac.es the 
asymptotic value of 0.5 for X/c * 2.25, indicating that the wakes from the 
adjacent blades are merged for X/c > 2.25. 

Figure 49 shows the influence of the rotor-stator spacing on the gust har- 
monic amplitude spectra. Tip vortex is not included in the computation of the 
spectra. The BPF level (m = l) is the highest for X/c — 0.5 and reduces mono— 
tonically as X/c increases. The BPF level for X/c = 5 is about 13dB lower com- 
pared to the BPF level for X/c =0.5. Since the wake width at X/c = 0.5 is 
the lowest, the fall-off rate of the spectra is the lowest for X/c = 0.5. As 
X/c increases, the wake width increases and the wake shape within a blade pas- 
sage resembles a sinusoidal function with a period equal to the blade-to-blade 
spacing. Hence, the contribution from the higher harmonics is less. There is 
a distinct change in the spectra as the semiwake width reaches the asymptotic 


104 



Normalized Semiwake Ncimalized Wake Centerline 

Width. 5/S Defect. W(j C /W 0 


• Rotor 55 

• 80% Design RPM 

• 10% Span from ""‘P 

• Cq = 0.1070 




figure 48. Axial Variation of Predicted Wake Centerline Defect and 
Semiwake Width. 


105 


• Rotor 55 

• 80% Design RPM 

• 10% Span from Tip 

• Neither Tip Nor Hub Vor'ex 

• Cp = 0.1070 



Harmonic Number, m 


Figure 49. Predicted Gust Harmonic Amplitude Spectra at 
Various Rotor/Stator Spacings. ? at 


value of 6/S = 0.5 for X/c > 2). This change in character is due to the merg- 
ing of wakes from the adjacent blades. Also note that once the wakes merge, 
the spectral shape does not change as X/c increases above 2. The slight 
reductions in the levels seen are due to a decrease in the centerline wake 
defect. Thus, increasing the rotor-stator spacing above a value at which the 
wakes from adjacent blades merge does not appreciably reduce the fluctuating 
lift component. This can be used as a design guide in the selection of the 
rotor-stator spacing, from the point of reducing the rotor wake interaction 
with the stator. 


5.2.2 Spanwise Variation in Rotor Wake Char ac ter is t ics 


This section deals with the influence of the spanwise variation of the 
geometric and aerodynamic parameters of the rotor on the rotor wake character- 
istics. Computations are done for Rotor 55 at 80% design rpm (see Reference 35 
for all the pertinent geometric and aerodynamic parameters). The gust har- 
monic amplitude spectra are computed with and without the tip vortex in order 
to study the influence of the tip vortex at different spanwise locations and 
downstream locations. 

Figure 50 shows the spanwise variation of rotor solidity (a r ), normalized 
axial spacing (X/c) and the rotor stagger angle ( Y r). The rotor-stator axial 
spacing shown in Figure 50 refers to the smallest spacing at which the rotor 
wake data were measured by Shaw and Balombin (Reference 6 ); it equals 0.54c at 
50% span. Figure 51 shows the spanwise variation of the work coefficient 
(Vq/Vt)> drag coefficient (C^), and rotor exit relative flow angle (g 2 ) at 
three rotor speeds (80%, 96%, and 115% design rpm). Note that only the drag 
coefficient near the tip shows some variation with the rotor speed. The drag 
coefficient is lowest for the 96% design rpra case and increases on either side 
of this rpm, vh^h is close to the design rpm. 

Figure 52 shows the predicted axial variation of the wake centerline 
defect (W dc /W 0 ), normalized streamwise distance (s/c) and the semiwake width 
( 5 /S), at three spanwise locations. The combination of the axial spacing 
(X/c) and rotor exit relative flow angle (& 2 ) at the three spanwise locations 
is such that the normalized streamwise distances at the three spanwise loca- 
tions are about the same. Thus, the wake centerline defect and the semiwake 
width are essentially dependent on the drag coefficient at these spanwise 
locations. Since the drag coefficient is lowest at the midspan location (50% 
span from tip), the normalized wake centerline defect and the semiwake width 
are minimum at this spanwise location. The wake centerline defects at 10% and 
90% span locations are about the same. The normalized semiwake width for the 
case of 90% span from tip is slightly higher than for the case of 10% span 
from tip. This is due to slightly higher streamwise distance and a slightly 
higher drag coefficient near the hub. 

Figure 53 shows the predicted gust harmonic spectra at the three spanwise 
stations (10%, 50%, and 90% span from tip) at three axial spacings without a 
tip vortex. The three axial spacings (X/c = 0.54, 1.23, and 1.77) refer to 
the rotor-stator spacings at raidspan at which rotor wake data were measured by 
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Shaw and Balombin (Reference 6). As X/c increases, the fall-off rate of the 
gust harmonic amplitude increases because of the growth of semiwake width. 

The fall-off rate for the 90% span from tip case is highest since its wake 
width is the highest. Note the reduction in BPF level as X/c increases at 
all three spanwise stations, because of the reduction in wake centerline 
defect. Also, at each axial location, the BPF level for the 50% span from tip 
location is lowest since the wake centerline defect is lowest for this span- 
wise location. 

Figure 54 shows the predicted gust harmonic spectra at the sa ae axial and 
spanwise locations as above, but including the tip vortex. The t agential 
location of the tip vortex is set at b t /S = 0.42, as in the case of data- 
theory comparison for Rotor 55 (see Subsection 5.1.1). Note there are dis- 
tinct changes in the spectral content for 10% span from the tip location (Fig- 
ures 53 and 54) at all three axial stations. The 2xBPF levels are higher than 
the BPF level for this spanwise location and the fall-off rate is reduced due 
to the contribution from the tip vortex. The spectra at 50% and 90% span from 
tip are not influenced by the tip vortex, since they are sufficiently far away 
from the tip vortex center. 

As shown in Subsection 3.3.1, the twist of the rotor blade from hub to 
tip results in a time delay (or phase lag) in the wake from hub portion imping- 
ing on the stator relative to the wake from the tip portion impinging on the 
stator. Figure 55 shows the predicted spanwise variation of the relative 
phase lag normalized by the angular spacing of the rotor for three rotor-stator 
spacings for Rotor 55 at 80% design rpra. Suppose that at t = 0, the wake from 
the tip region of a particular rotor blade m A" impinges on the stator. For the 
case of X/c =1.77 (see Figure 55), by the time the wake from the hub region 
(95% span from tip) from the same rotor blade n A M impinges on the stator, 
approximately 4.6 rotor blade passages would have gone past the stator. Thus, 
in Figure 55 we see that as X/c increases, the phase lag increases. The extent 
of spanwise variation of the phase lag depends on the amount of twist in the 
rotor blade: the larger the twist, the more the phase lag. 


5.2.3 Influence of Rotor Solidity 

Rotor 55 at 80% design rpra is again chosen to study the influence of rotor 
solidity. The streamline at 10% span frou tip is considered so that influence 
of rotor solidity on the stator upwash gust harmonic spectra can be studied 
with and without tip vortex. The solidity for Rotor 55 at 10% span from tip 
equals 0.905. The range of variation in the solidity studied is ±33% of the 
nominal valve of 0.905. The rotor solidity is varied by changing the rotor 
aerodynamic chord while keeping the same blade-to-b 1 ade spacing. The wake/ 
vortex ccnputations are performed at the same axial distance from the rotor 
trailing edge for these parametric studies. The X-locations correspond to the 
locations where Shaw and Balombin measured the wake data (Reference 6). As 
solidity increases, the rotor aerodynamic chord increases. As the rotor aero- 
dynamic chord increases, X/c decreases for constant X values. 
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Figure 56 shows the influence of rotor solidity on the axial variation 
of the wake centerline defect and the semiwake width. As the rotor solidity 
increases, the wake centerline defect increases, and the semiwake width 
decreases at the same X-location due to reduced X/c value. Thus, the semi- 
wake width reaches the limiting value of 0.5? for lower solidity case at a 
smaller axial distance from the rotor trailing edge, compared to the higher 
solidity case. 

Figure 57 shows the predicted influence of rotor solidity on the gust 
harmonic spectra at three axial locations from the rotor trailing edge without 
the tip vortex. Since the case of the highest rotor solidity has the highest 
value of wake centerline defect at X/r ann = 0.18, the BPF level for it is 
highest. Also, since the semiwake width is lowest for this case, the fall-off 
rate is also the lowest of the three solidity cases. The same trend is noticed 
at X/r ann * 0.71. For X/r ann = 1.43, the semiwake width has reached its limit- 
ing value of 0.5S for all three rotor solidities. Thus, the gust harmonic 
spectra for all three rotor solidities are almost identical. The slight dif- 
ferences are due to differences in the wake centerline defect for them at this 
X-location . 

Figure 58 shows tne predicted influence of rotor solidity on the gust 
harmonic spectra at these same three axi il locations, but this time including 
the tip vortex. The tangential location of the tip vortex for this study is 
set at b t /S - 0.5. Note that the significantly different character in the 
spectra for this streamline is due to the tip vortex. At X/R anr = 0.18, the 
2xBPF and 4xBPF levels are significantly higher than the BPF and 3xBPF levels, 
respectively. Since the tip vortex is exactly in the raidpassage, the stream- 
wise defects due to the wake and the tip vortex (in the relative frame of 
reference) yield two well-defined defects within one blade-to-b lade spacing. 
This increases the contribution of even harmonics of BPF and lowers the con- 
tribution of odd harmonics of BPF. The tip vortex dominates the spectra at 
this spanwise location. A clear trend with rotor solidity cannot be extracted 
in the presence of tip vortex at streamlines close to the tip. As X/r ann 
increases, the tip vortex decays and the relative dominance of 4xBPF over 
3xBPF is not seen. However, 2xBPF levels are still higher than BPF levels, 
even at the X/r ann = 1.43. 


5.2.4 Influence of Section Drag Coefficient 


A parametric variation of the section drag coefficient (C D ) and its 
influence on the wake centerline defect, semiwake width, and the stator upwash 
gust harmonic spectra are discussed in this subsection. As in Subsection 
5.2.3, the streamline at 10% span from tip for Rotor 55 at 80% design rpm is 
considered. The nominal value of Cq computed from the measured profile loss 
coefficient equals 0.1070 ,Cd was varied ±50% from this nominal value and 
its influence was studied. 

Figure 59 shows the axial variation of the predicted wake centerline 
defect and semiwake width for the three values of C D . The influence of C D 
on the wake centerline defect and semiwake width, in particular, is seen to 


115 







Figure 58. Predicted Influence of Kotor Solidity on Gust Harmonic Spectra at Three Axial 
Locations f rom Kotor Trailing Edge (Tip Vortex Included). 
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' * 2 ‘ ^ .Influence ot R otor Loading at Constant Tip Speed 

The pressure ratio of a rotor is the indicator of rotor loading. For a 
single stage rotor, such as Rotor 55, the rotor loading does not va^y notice- 
ably along a constant speed line (see Reference 35). This is not the case for 
rh l T *a rotors . To study the effect of rotor loading at constant tip speed 
the second stage of Rotor 67 at 1002 design -pm was chosen (see Reference 37) 
mce it exhibited significant variation in rotor pressure ratio. 

ar rh I 1 ? Ure 61 shows the variation of the work coefficient, lift coefficient 
at the tip and section drag coefficient with the rotor total pressure ratio 

°l ll*** 2 H f R °T 67 - at 102 3Ran fr3 “ tLp ’ at 100% desi S n rpm. The needed 
geometric and aerodynamic data were taken from Reference 37. The tip lift 

coefficient and the work coefficient are seen to increase rapidly with the 

rotor total pressure ratio. The drag coefficient is fairly Hat <^e would 

anticipate an increase in drag coefficient with increase in rotor total pr^- 

Roto/bWR f exam ^ atl °" of th * performance map of the second stage of 
Rotor 67 (Reference 37) indicates that the adiabatic efficiency of the rotor 
decreases rather sharply with decrease in rotor total pressure ratio This 
has resulted m a fairly flat C D variation with the rotor total pressure 
ratio. Figure 62 shows the variation of the section drag coefficient (again) 
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, ,, / e wa e centerline defect and the semiwake width are computed at the 
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s little variation with the rotor total pressure ratio, the normalized 

t^rolo^ 1 ! 11 ? defeCt and the semiwake width also show little variation with 
lolorlol*? pressure ratio. Figure 63 shows the predicted influence of 

: tal ,r SS T ratl ° at 1002 design rpm 0n che 8 ust h a ™°"ic spectra, 

with and without the tip vortex. Because the wake centerline defected the 

semiwake width are almost constant over the range of the rotor total pressure 
ratio considered, the gust harmonic spectra for the case of no tip voltex show 
negligible variation with the rotor total pressure ratio. However the situa- 
tion is quite different when tip vortex is included in the compulsion of gult 

l ™8?\TTlnn ^ , CaSe ° f thG high6St r0t ° r t0tal ratio (?" l 

the 2*BPP level TTT h tip vortex is noticeably strong in the sense that 

the 2xBPF level is higher than the BPF level. As the rotor total pressure 

10 decreases, the influence of tip vortex is observed to decrease. The 
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gust harmonic spectrum for P r = 1-281 with and without the tip vortex is not 
significantly different. The reason for the above behavior is that the lilt 
coefficient at the tip increases with the rotor loading. Therefore, for the 
same fraction of lift lost to the tip vortex, (k) ti]? , the circulation 

of the tip vortex increases (see Subsection 3.2.1) and the influence of tip 
vortex on the stator upwash velocity and the corresponding gust harmonic 
spec trum increases. 


5.2.6 Influence of Inviscid Velocity Gradient 

The influence of the wake asymmetry in terms of a prescribed velocity 
gradient on the rotoi wake and vortex characteristics was studied m this 
section. For the sign convention adopted for rotor wakes (see Figure 1), 
a negative value for the inviscid velocity gradient typically represents the 
fan rotor wakes (see Figure 15) because in such a case the relative velocities 
on the pressure side are lower than on the suction side. 

The streamline at 10% span from the tip of Rotor 55 at 80% design rpm was 
chosen for this study. Figure 64 shows the predicted influence of three values 
for the inviscid velocity gradient (WTIV = -1.0, 0.0 and +1.0) on the tangen- 
tial distribution of the relative velocity normalized by wheel speed, over one 
blade passage, at three values of rotor-stator spacings. The tip vortex is 
not included in these predictions. The extent of asymmetry decreases as X/c 
increases. Note the excessive asymmetry in the relative velocity at X/c 0.5. 
The tangential profile for the relative velocity for WTIV = -1.0 is a mirror 
image of the profile for WTIV = +1.0 for all X/c locations about the wake 

centerline . 

Figure 65 shows the influence of inviscid velocity gradient (WTIV = -1.0, 
0.0 and +1.0) on the tangential distribution of the normal perturbation veloc- 
ity (that is, gust velocity) w ch and without the tip vortex at X/c - 0.5- The 
tangential distribution of the gust velocity is shown for the case of a tip 
vortex at b t /S *0.5 and 0.42, as well as for the case of no tip vortex. 

In the case of no tip vortex, the tangential profile of the gust velocity for 
WTIV * -1.0 is a mirror image of the profile for WTIV = +1.0 about the wake 
centerline. The tip vortex modifies considerably the tangential distribution 
of the gust velocity. For the case of the tip vortex located at b t /S - 0.5, 
the tangential profile of the gust velocity for WTIV = -1.0 is still seen to 
be an almost mirror image of the profile for WTIV = +1.0. This happens because 
the tip vortex is exactly at the midpassage. However, when the tip vortex is 
located at b t /S - -0.42, the tangential profile of the gust velocity for 
WTIV * -1.0 no longer resembles the profile for WTIV = +1.0. Also, there was 
a dip in the gust velocity at the tangential location of the tip vortex for 
all three values of WTIV. Thus, it is seen that the inviscid velocity gradient 
in combination with the tip vortex can yield complex tangential profiles of 
the gust velocity. 

Figures 66 through 68 show the influence of the inviscid velocity gradi- 
ents and tip vortex on the gust harmonic amplitude spectra at three values of 
rotor/stator spacing. Figure 66 shows the influence of the inviscid velocity 
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gradient on the gust harmonic amplitude spectra at three values of rotor/stator 
spacing for the case of no tip vortex. Since the tangential profiles of the 
gust velocity for WTIV * +1.0 and -1.0 are mirror images about the wake center- 
line, their Fourier decomposition yields the same Fourier coefficients and 
identical spectra. For the case of WTIV = 0.0 (no asymmetry), the higher har- 
monics fall off more rapidly than in the case of asymmetry. Figure 67 shows 
the influence of the inviscid velocity gradient on the gust harmonic amplitude 
spectra at three values of rotor/stator spacing for the case of tip vortex 
located at b t /S * 0.5. Since the tangential profiles of the gust velocity 
for WTIV * +1.0 and -1.0 are almost mirror images of each other, the spectra 
for WTIV - +1.0 and WTIV = -1.0 are identical. We also notice significant 
dominance of the tip vortex at X/c = 0.5 for all three values of WTIV. Figure 
68 shows the influence of the inviscid velocity gradient on the gust harmonic 
amplitude spectra at three values of rotor-stator spacing for the case of tip 
vortex located at b^/S - 0.42. The inviscid velocity gradient influences 
the BPF levels more than the harmonics of BPF at all three rotor/s t ncor spac- 
ings. At X/c = 0.5 and 2.0, the 2xBPF levels are higher than the BPF levels 
for WTIV = 0.0 and +1.0, whereas for WTIV = -1.0, ^he BPF level is slightly 
higher than the 2xBPF level. This indicates that the negative inviscid veloc- 
ity gradient (WTIV * -1.0) can offset the influence of tip vortex by pushing 
the 2xBPF level below the BPF level. Thus, a combination of wake asymmetry 
and tip vortex can yield complex shapes for the gust velocity spectra. 


5 • 2 • 7 Influence of Tip Clearance and Tangential Location of Tip Vortex 

It has been shown in the previous subsections that the tip vortex can 
have a significant effect or the stator upwash velocity spectrum. This sub- 
section deals with the influence of two parameters of the tip vortex (tip 
clearance and tangential location of tip vortex) on certain parameters of the 
vortex structure on the stator upwash gust harmonic spectra. The stream- 

line located a } , 10% span from tip for Rotor 55 at 80% design rpm is chosen for 
this parametric study. 


69 shows the predicted variation of normalized radius of vortex 
core, streamwise velocity defect of vortex core, and circulation of vortex with 
rotor-stator spacing for five values of tip clearances. The mathematical 
modeling of these parameters is described in Subsection 3.2.1. The nominal 
value representing the tip clearance of Rotor 35 is / ( c ) t = 0.009. Tip 
clearance has been varied around this nominal value. The tip vortex radius 
increases with tip clearance, as well as with rotor-stator spacing. The vortex 
radius is modeled to increase as (s/c + The streamwise velocity 

deficit of the vortex core increases with the radius of vortex core, which in 
turn increases with tip clearance. The streamwise velocity deficit of the 
vortex core is modeled to decay as l/(s/c+l) (see Subsection 3.2.1). Figure 69 
also shows the decay of the circulation of the tip vortex with axial distance 
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Previous sections have demonstrated that the presence of a tip vortex at 
the mid blade passage yields anocher defect in the tangential profile of rela- 
tive velocity in addition to that of the rotor wake within one blade passage. 
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Hence, 2xBPF levels of the gust velocity harmonic spectra have tended to be 
than the BPF levels. A more detailed study of this phenomenon is 
reviewed here. Figure 70 shows the predicted influence of tip clearance on 
the gust harmonic amplitude of BPF and 2xBPF tones at three rotor-stator 
spacings for the streamline at 10% span from the tip for Rotor 55 at 80% 
design rpm. The tangential location of the tip vortex for these computations 
is set at b t /S = 0.5. The 2xBPF levels are larger than the BPF levels at 
all three rotor-stator spacings. For X/c = 2 and 4, the ^dB between 2xBPF 
and BPF levels decreases, indicating decay of the vortex at these distances. 
Also, we notice sharp reduction in BPF level in the vicinity nf T /(c)^£ = 
0.012, at X/c = 0.5. For values of t/(c)^£p on either side, the &dB between 
2xBPF and BPF levels decreases. This behavior can be explained by examining 
the tangential distribution of the gust upwash velocity (Figure 71). Figure 
71 shows the tangential distribution of the gust upwash velocity over two 
blade passages for r/(c) ci = 0.003, 0.012, and 0.02, at X/c = 0.5. As the 
tip clearance increases, the peak defect due to the tip vortex (that is, 
defect at x/S = 0.5 or 1.5) increases compared to the rotor wake (that is, 
defect at x/S — 0.0 or 1.0 or 2.0). At t/(c)^£p = 0.012, the peak defect 
created Dy the tip vortex is almost identical to that due to the rotor wake. 
Thus, for t/(c) t £p = 0.012, the gust upwash velocity tangential distribution 
shows two almost identical waveforms within one blade passage. Hence, the 
F. •'•ier analysis of such a waveform shows a negligible contribution from the 
BPF hacn.-mic amplitude relative to the 2xBPF amplitude. On either side of 
this tip clearance, the peak defects due to tip vortex and rotor wake are not 
identical, even though there are two waveforms within one blade passage; thus, 
the AdB between 2xBPF and BPF levels decreases. This is the reason why the 
BPF level tends to approach 2xBPF tone level for increasing tip clearances at 
other axial locations. 

Figure 72 shows the predicted influence of tip clearance on the gust har- 
monic amplitude of BPF and 2xBPF at three rotor-stator spacings for the case 
of b[/S = 0.45, for the streamline at 10% span from tip for Rotor 55 at 80% 
design rpm. Since the tip vortex is not exactly at the middle of the blade 
passage, there is not a sharp fall in the BPF level at X/c =0.5, as was 
noticed in Figure 70, where the tip vortex was exactly in the middle of the 
blade passage. Figure 73 compares the predicted tangential distribution of 
the gust upwash velocity for b t /S = 0.5 and 0.45, for T /(c) t £ p = 0.012. For 
T / (c)tip ~ 0*012 and b^/S = 0.5, almost two identical waveforms exist within 
one blade passage. When the tip vortex is at b t /S = 0.45, there are no two 
identical waveforms within one passage, even though the peak defects due to 
rotor wake and tip vortex are almost the same. The relative variation between 
2xBPF and BPF ampliLuc.es at X/c = 2 and 4 as a function of tip clearance for 
b t /S = 0.45 follows that of b t /S =0.5 (Figures 70 and 72). Thus, one notes 
that the behavior of the upwash velocity spectrum near the tip region is very 
sensitive to the tip clearance and tangential location of the tip vortex. 

Figure 74 shows the predicted influence of the tangential location of tip 
vortex on the gust upwash harmonic amplitude spectra at three rotor-stator 
spacings for the streamline at 10% span from tip for Rotor 55 at 80% design 
rpm. The gust upwash harmonic amplitude spectrum for the case of no tip vor- 
tex is also included for the sake of reference. At all three rotor-stator 


134 


Gust Upwash Velocity/Wheel Speed 


• Rotor 55 

• 80% Design RPM 

• 10% Span from Tip 


• bt/S = 0.5 

• X/c = 0.5 



Tangential Spacing, x/S 


Figure 71. Predicted Tangential Distribution of the Gust 
Upwash Velocity for Three Values of Tip 
Clearance . 


135 





Gust Up wash Velocity /Wheel Speed 


Rotor 55 

80% Design RPM 

10% Span from Tip 




Tangential Spacing, x/S 


Figure 73. Predicted Tangential Distribution of the Gust 

Upwash Velocity for Two Tangential Locations 
of the Tip Vortex. 




8P ‘ (paadg jaaijM aATtjpja^) 
apn^x^duiv otuouuph xsnn 






M 


spacings , the 2xBPF levels are higher than the BPF level when the tip vortex 
is located at b c /S * 0.5. In fact, at X/c = 0.5, the 4xBPF level is higher 
than the 3xBPF level, for b c /S =0.5. As the tip vortex is moved awav from 
the blade midpassage to b t /S = 0.25, the 2xBPF level is no longer higher. 

In fact, the odd multiples of BPF are higher than the subsequent even multi- 
ples of BPF. When the tip vortex is located right at the wake centerline of 
b t /S - 0.0, the fall-off rate of the harmonics is very slow. For the case of 
no tip vortex, the spectrum shows monotonic fall-off with harmonic number. 

^ ~ ^ an< * the fall off rate is quite high, indicating sufficiently 

large semiwake widths, so that the wake profile within one blade passage resem- 
bles a sinusoidal wave with blade passage as its period. Also, at large spac- 
ings (X/c - 2 and 4), there is a significant difference between the no tip vor- 
tex case and cases with a tip vortex at any tangential spacing. This implies 
uhat, for streamlines close to the tip, no matter what the tangential location 
of a tip vortex, it can significantly modify the gust upwash velocity spectrum. 


Figure 75 shows the predicted influence of the tangential location of the 
tip vortex on BPF and 2xBPF levels at three rotor-stator spacings in order to 
highlight some of the features presented in Figure 74. We note that at X/c = 
0.5, the 2xBPF level is higher than the BPF if 0.36 S b t /S £ 0.5. At X/c « 

2.0, the 2xBPF level is higher than BPF if 0.41 £ b t /S £ 0.5. For X/c = 4.C 
the 2xBPF level is higher than BPF if 0.44 S b t /S £ 0.5. This indicates that 
as k/c increases, the existence of two similar waveforms within one blade pas- 
sage (one due to rotor wake and another due to tip vortex) is restricted to 
narrower regions about the blade midpassage where the tip vortex ought to 
exist. This happens because the tip vortex and the rotor wake decay at dif- 
ferent rates. At X/c = 0.5, there is a sharp drop in the 2xBPF level when 
b t /S = 0.25. As shown in Figure 74, odd multiples of BPF are predominant at 
this location. An explanation for this and other features observed in Figure 
74 can be found by examining the tangential distribution of the gust upwash 
velocity (see Figure 76). Figure 76 shows the predicted tangential distribu- 
tion of the gust upwash velocity for three tangential locations of the tip 
vortex (b t /S — 0.5, 0.25, and 0.0) and for the case of no tip vortex, for 
X/c — 0.5. For b t /S — 0.5, there are almost two identical waveforms within 
one passage. Therefore, even multiples of BPF show higher levels than the pre- 
ceding odd multiples of BPF (see Figures 74 and 75). For b t /S * 0.25, the 
defects due to tip vortex and rotor wake interact, and the resulting waveform 
is not symmetric as it was for b t /S = 0.5. This causes the odd multiples of 
BPF to increase more than succeeding even multiples of BPF. When b t /S = 0.0, 
the defects due to wake and vortex yield a sharp upwash velocity profile at the 
blade centerline, similar to a step function. The Fourier transformation of an 
impulse function is a constant (that is, a step function can be constructed by 
summing infinite harmonics of constant amplitude). The gust upwash harmonic 
amplitude spectrum for bj./S = 0.0 indeed shows the lowest fall-off rate with 

V» a me r*» i • *nK T« fUe h,’_ t , , 
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profile shows one defect due to the wake, and thus the higher harmonics fall 
mono t onical ly with the harmonic number. 
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5-2.8 Influence of Velocity and Length Scales on Turbulence Spectrum 


This subsection deals with the influence of axial and transverse velocity 
and length scales on the axisymmetr ic turbulent velocity spectrum model, which 
was described briefly in Subsection 3.3.2. 

The streamline at 10% span from the tip of Rotor 55, operating at 80% 
design rpm, was chose i for this parametric study. The predicted axial and 
tangential turbulence velocities at X/c = 0.54 at raidspan equal 0.0333 x 
v wheel> anc * the pcedictec radial turbulence velocity equals 0.0394 x V w ^ ee ^ . 
Since there is not enough information available on axial and transverse 
length scales, they are assumed to be equal, and are prescribed as fractions 
and multiples of the semiwake width. The above velocity and length scales 
yield a value of 1.8 for the turbulence asymmetry parameters, e (see Sub- 
section 3.3.2). 

Figure 77 shows the predicted influence of varying the axial and trans- 
verse velocity scales while holding the axial and transverse length scale 
constant at L a /5 = L t /<5 =0.5, and also the predicted influence of varying 
the axial and transverse length scales while holding the axial and transverse 
velocities constant at the nominal values of u f /V^ ee ^ = ^-0333 an< * 

w '/ v wheel ® 0.0394. The bandwidth for these predictions equals 0.5 BPF . 

The axial and transverse velocity and length scales are varied so that e is 
coustaut at a value of 1.8. We notice that doubling the velocity scales 
increases the broadband levels by about 6 dB and the shape of the broadband 
spectrum is not altered. However, varying the length scale alters the levels 
and the shape of the broadband turbulence spectrum. At smaller length scales, 
the broadband levels fall more slowly with increasing frequency. 
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6.0 CONCLUSIONS AND RECOMMENDATIONS 


The principal objective of this research program was to establish a veri- 
fied rotor wake and vortex model for specific application to fan and compres- 
sor rotor-stator interaction and resulting noise generation. The program 
called for the development of an empirical model for rotor wake and tip vortex 
characteristics, along with a turbulence spectrum model to predict the broad- 
band spectrum of the stator gust upwash. 


6.1 CONCLUSIONS 

Certain empirical rotor wake and turbulence relationships were developed 
using the low speed rotor wake data of the Pennsylvania State University 
rotors . A tip vortex model was developed by replacing the annulus wall with a 
row of image vortices. An axisymmetric turbulence spectrum model, developed 
in the context of rotor inflow turbulence, was adapted for predicting the tur- 
bulence spectrum of the stator gust upwash. An interesting conclusion, arrived 
at while correlating the tangential distribution of turbulent velocities, is 
that the tangential profile of turbulent velocities is fairly flat, whereas 
the tangential profile of the turbulent intensity is Gaussian. 

The principal conclusions drawn from the data-tneory comparisons performed 
in this research program are as follows: 

• The agreement between the predicted and measured gust harmonic ampli- 
tude spectra for the JT15D fan rotor at various span locations for 
four rpm's along the engine operating line is excellent. 

• The empirical rotor wake and vortex model has been shown to predict 
very well the observed trends both in the gust harmonic spectra and 
the turbulence spectra for Rotor 55 at three speeds. 

• The predicted and measured wake centerline defect and semiwake width 
for Rotor 67 are in reasonable agreement for axial distance greater 
than about 0.4 chord, out do not agree close to the rotor trailing 
edge. The reason for the disagreement is believed to be the shock 
boundary layer interaction over the rotor blades of this high speed 
rotor, which significantly modifies the wake flow close to the trail- 
ing edge. The prediction model has been developed using low speed 
rotors which do not exhibit the shock boundary layer interaction. 


The table below gives the range of certain key parameters of the 
df.a base employed to develop the empirical relations for the pre- 
diction model and the range over which the prediction model has been 
verified to yield reasonable data-theory comparisons. 
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Parameter 


Axial Spacing at 
Midspan (X/c) 

Tip Speed 

Rotor Solidities 


Range of Parameter in 

the Data Base 

0.01 to 1.664 
95 fps to 168 fps 
0.68 to 1.28 


Range of Parameter 
Over Which the 
Model is Verified 

0.4 to 1.8 

560 fps to 1000 fps 

0.9 to 1.772 


A significant portion of this effort was devoted to parametric studies of 
the verified rotor wake and vortex prediction model. Some of the interesting 
conclusions of these parametric studies are given below: 

• Increasing the rotor-stator spacing above a value at which the merg- 
ing of wakes from the adjacent rotor blades takes place does not 
yield any appreciable reduction in the fluctuating lift component. 
This can be used as a desigr guide in the selection of rotor-stator 
spacing from the point of view of reducing the rotor wake interac- 
tion with the stator. 


• Tip vortex yields a streamwise defect just like the wake from the 
rotor blade, and also induces a significant amount of azimuthal and 
spanwise velocity components. By virtue of this, the tip vortex can 
significantly influence the gust spectrum, particularly near the tip 
region. The presence of a tip vortex typically results in higher 
levels for 2*BPF compared to BPF , which is essentially due to the 
presence of two defects within one blade passage. 

• The tip clearance is a critical parameter in determining the strength 
of the tip vortex, in that the higher the tip clearance, the higher 
is the fraction of lift lost to the tip vortex. At a certain value 
of the tip clearance, the streamwise defect caused by the tip vortex 
equals that caused by the rotor wake. If the tip vortex were to be 
located at exactly the midpassage in such a case, the rotor exit 
flow has two identical defects within one blade passage, which 
results in a gust harmonic spectra dominated by even harmonics of 
BPF. 

• The tangential location of the tip vortex has an important bearing 
on the gust harmonic spectra. If the tip vortex is exactly in the 
blade midpassage, even harmonics of BPF dominate the gust spectra; 
if the tip vortex is exactly at 1/4 blade passage, odd harmonics of 
BPF dominate. If the tip vortex is right at the wake centerline, 

the combination of the wake and vortex defects yields a narrow defect 
which results in increased contribution to the gust harmonic spectra 
by all harmonics of the BPF. 


* 
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• Rotor loading (rotor pressure ratio) causes asymmetry of rotor wakes 
on the pressure and suction sides of the rotor blades. A simple 
model of the rotor wake asymmetry in terms of a prescribed inviscid 
velocity gradient was developed to study this aspect. Wake asymmetry 
by itself did not have much effect on the gust harmonic spectra even 
at smalL rotor/stator spacings. However, a combination of asymmetric 
wake and tip vortex results in complex gust harmonic spectra. 

• The section drag coeff ci’nt is normally minimum at the midspan loca- 
tion of a rotor blade, s, :he wake defect and the semiwake width are 
lowest in the midspan location. 

• The section drag coefficient is normally minimum at the design rpm, 

so the wake defect and the semiwake width are lowest at the design 
rpm. 

• Higher rotor solidity (obtained by increasing the rotor chord and 
keeping the same blade-to-blade spacing), results in larger wake 
defects and smaller semiwake widths at a prescribed distance from 
the rotor trailing edge. Hence, the stator upwash spectra (without 
tip vortex) shows higher absolute level and lower fall-off rate with 
the harmonic number with higher solidity, as long as the wakes from 
the adjacent blades have not merged. Once the wakes from the adja- 
cent blades merge, the rotor solidity has very little influence on 
the stator upwash spectral fall-off rate with the harmonic number. 

A ciear trend with rotor solidity cannot be extracted in the pres- 
ence of a tip vortex at streamlines close to the tip because the 
tip vortex dominates the stator upwash spectra at these streamlines. 

In summary, this research program has established and verified a rotor 
wake and vortex model. The parametric studies have given a valuable insight 
into the influence of various geometric and aerodynamic parameters on the 
stator upwash gust harmonic spectra. The computer program requires about 10 
seconds of CPU time per case on the GE Honeywell-6000 system and thus is 
ideally suited for extensive parametric and preliminary design studies. 


6.2 RECOMMENDATIONS 


There are several possible improvements for the above developed rotor wake 

and vortex model. Some of the salient improvements recommended are discussed 
as follows. 

• In recent times, the counterrotating propeller concept has been 

shown to be a very attractive way to reduce specific fuel consump- 
tion. It is anticipated that there will be strong rotor-to— rotor 
interactions because of the small axial spacings between the two 
counterrotating propellers, possibly introducing high noise levels 
and/or forced vibrations. The rotor-stator wake and tip vortex 
interaction model developed here could be generalized to the case of 
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the counterrotating propellers to aid in the understanding and pre- 
diction of the rotor-rotor interaction noise mechanisms. It is rec- 
ommended that such a study be undertaken. 

The tangential location of the tip vortex could not be ninned down 
by correlating the existing experimental data of the wave forms of 
the upwash component in the tip region, and hence it is prescribed 
as an input parameter. At present this parameter is somewhat arbi- 
trary in the sense that we have to adjust this parameter to get 
reasonable agreement between data and theory. The detailed paramet- 
ric studies show the sensitivity of the gust harmonic spectra to 
this parameter. A detailed experimental survey of the tip region 
waveforms of the upwash component at several axial stations for 
various rotors at different operating speeds is recommended. A cor- 
relation for predicting the tangential location of the tip vortex in 
terms of tip clearance, rotor speed, rotor loading at the tip, etc., 
can then be developed. 

The rotor wake momentum integral analysis (Section 4.2) offers the 
ability to predict rotor wake structure, accounting for many more 
influential parameters with much less empiricism than the current 
correlation approach (Section 3.1). The choice of blade spacing, 
instead of rotor aerodynamic chord, as the normalizing dimension for 
the semiwake width, is reasonable for the high solidity rotors being 
considered in the current correlation approach. However, for low 
solidity rotors, rotor aerodynamic chord presumably is the correct 
normalizing dimension. The rotor wake momentum integral analysis 
will clarify this issue and will aid in the proper choice of the 
characteristic dimension for rotor wakes of widely varying solidities. 

As the parametric studies have indicated, the stator upwash harmonic 
spectrum is influenced significantly by the combined presence of the 
tip vortex and the inviscid velocity gradient (that is, rotor load- 
ing). At the present time, the inviscid velocity gradient is a pre- 
scribed input and is not correlated with rotor loading. A cascade 
flow analysis is recommended so that the inviscid velocity gradient 
can be correlated with the rotor loading. 

Though the prediction model can handle both tip and hub vortices, 
the rotor data surveyed did not give any quantitative information 
useful in developing empirical relations to compute the hub vortex 
strength and radius. It is recommended that, if and when such data 
become available, the empirical relations for hub vortex strength 
and radius be finalized. 
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8.0 NOMENCLATURE 


a 

AR 

B 

b 


BPF 


C L 


DoD 

DTIC 

GE 

h 


IGV 

k 

5 (k*, kp, kg ) 

L 

l 

La 

L d, L s 


Radius of the vortex core 
Aspect ratio (span/chord) 

Number of rotor blades 

Tangential distance of the center of the vortex core from 
the wake centerline (see Figure 26) 

Blade Passing Frequency 

Rotor aerodynamic chord 

Ambient speed of sound 

Drag coefficient 

Lift coefficient 

m-th Complex Fourier coefficient of the gust upwash 
ve loc ity component 

Department of Defense 

Defense Technical Information Center 

General Electric Company, USA 

Rotor blade span 

/=! 

Inlet Guide Vane 

Fraction of suction lift lost to the vortex 
Wave number vector 
Turbulence length scale 
Characteristic length scale 
Axial turbulence length scale 

Semiwake width at half the depth on the pressure and 
suction sides of the wake, respectively 
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L t 


m 




M,N, 


n 


NACA 

NASA 

0 (#) 

P 

P r 

R 

r 

R SJ RN 
r * 

R/S 

S 

s 

U 

(u\ V 1 , w 1 ) 

u c 

u d 

u dc 


Transverse turbulence length scale 
m = 1 3PF 


Harmonic number 


m = 2 2 x BPF , etc. 


Axial Flow Mach No- 

Parameters associated with vortex (see Equations 27 through 
29 and Figure 26) 

Normal to the streamwise coordinate (in the tangential 
direction) 

National Advisory Committee for Aeronautics 
National Aeronautics and Space Administration 
Order of magnitude 
Static Pressure 
Pressure ratio 

Domain of the vortex core (see Figure 26) 

Radial distance 
Reynolds Number 

Distance from the vortex center 
Rotor-s tator 
Blade-to-b lade spacing 

Streamwise distance from the rotor trailing edge 
Axial velocity (see Figure 30) 

Axial, tangential, and radial component of turbulent (rms) 
ve loci ty 

Velocity defect at centerline 

Axial component of wake velocity defect 

Axial component of the wake centerline velocity defect 


u 


2 


o 

u 4 - 

t 


Axial and transverse components of turbulence intensity 


152 


iiiiifeii 


V 


Noruial velocity component (see Figure 31) 

Velocity vector in the stationary frame of reference 


r 

! 

r 


v ref 

v t 

^wheel 


W 

W 

«dv 

w d (y') 
w m 

Wi 


w r I 

( ) 


Perturbation velocity vector in the stationary frame of 
ref erenc e 

Upwash component of the perturbation velocity 

Streaarwise component of the perturbation velocity 

Reference velocity used in evaluating dB level for upwash 
velocity harmonic content 

Tip speed 

Wheel speed at a given radial location 
Tangential velocity (see Figure 30) 

Tangential component of the flow velocity at rotor exit 
(assuring no swirl at rotor inlet) 

Velocity vector in the rotating frame of reference 

Total velocity in the rotating frame of reference 

Maximum defect in the strearawise velocity of the vortex 
core in the rotating frame of reference 

Velocity defect in the streamwise component of velocity at 
a distance y* from the vortex center in the rotating frame 
of reference 

Mean fl-w velocity through the blade passage in the rotat- 
ing frame of reference 

Main stream flow velocity at the rotor inlet in the rotating 
frame of reference 

Work coefficient 

Tangential velocity induced by the ( ) vortex in the 

rotating frame of reference 

Radial velocity induced by the ( ) vortex in the rotating 

frame of reference 
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* 


f 


WTIV Inviscid velocity gradient normalized by wheel speed 

X Axial distance from the rotor trailing edge 

x > y Coordinates of the unwrapped annulus (see Figure 26) 

y) Cascade plane (see Figure 30, Section 3,3.2) 

(x 0 , Xo> Fluid fixed coordinates (see Section 3.3.2) 

AX Axial distance between rotor trailing edge and 1/4 chord 

point of stator as a point at which the gust velocity is 
being evaluated 


GREEK SYMBOLS 


a 


B 

81 

82 

r 

YR 

YS 

6 

«ij 


n 

e 

Xr * X s 

u 

v 


V 


T 


Flow angle of the absolute velocity vector to the axial 
direction 

Flow angle of the rotor relative velocity vector to the 
axial direction 

Rotor relative inlet flow angle 

Rotor relative exit flow angle 

Circulation 

Rotor stagger angle 

Stator stagger angle 

Semiwake width 

Kronecker's delta, 5 — = 1 if i = j 

* 0 if i + j 

Normalized tangential distance [ T ) * x/(6/2)] 

Momentum thickness in a two-dimensional boundary layer 
Rotor or stator stagger angles, used in Section 4.2 
Turbulence velocity scale (see Equation 104) 

Kinematic viscosity (laminar) 

Eddy viscosity 
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v t 

p 


T T 

rn> sn 


x, o 


(O 

x c 



A<|> 

X 

a 

CJ 

u 


Turbulent equivalent kinematic viscosity 
Density (fluid density) 

Rotor solidity (rotor aerodynamic chord/b lade-t o~b 1 ade 
spacing) 

Tip clearance 

Shear stresses in radial and streamwise directions 

Turbulent shear stress 

Axial component of turbulent intensity 

Axial component of turbulent intensity in the free stream 

Axial component of turbulent intensity at the wake center- 
line 

Angie from the rotor blade stacking axis to the 1/4 chord 
poin 4- of stator or a point at which the gust velocity is 
being evaluated 

Phase of the complex nr-th order Fouri^~ coefficient, c m 
Turbulence upwash velocity spectrum 

Angle from the rotor blade stacking axis to the trailing 
edge of the rotor 

♦ “ *o 

Turbulence wavenumber (see Equation 79) 

Angular velocity of the rotor 
Angular velocity of the vortex 
Turbulence spectrum frequency 


Subscripts 

a Refers to axi al/ tangential turbulent length scales 

ann Refers to annulus wall 

c Refers to centerline 

d Refers to wake defect 
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dc 


h , hub 

^ hub^ vtx 
n 


o 

ps 

s 

ss 

t 

t, tip 
( ti P ) vtx 


Refers to defect at the wake centerline 
Refers to edge 
Refers to hub 

Refers to hub vortex (particularly for circulation) 

Refers to normal to the streamwise direction (in the 
tangential direction) 

Refers to free stream condition 

Refers to pressure side of rotor blade 

Refers to the streamwise direction 

Refers to suction side of rotor blade 

Refers to transverse (radial) turbulent length scale 

Refers to tip 

Refers to tip vortex (particularly for circulation) 
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